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I .   INTRODUCTION 

A.   THESIS  SUMMARY 

The  objective  of  this  thesis  is  to  develop  a  software 
package  to  beamform  acoustic  signals  used  in  ocean  acoustic 
tomography.  The  goal  of  tomography  signal  processing  is  the 
precise  measurement  of  acoustic  travel  time.  Sound  speed  is 
a  well  understood  function  of  temperature,  pressure  and 
salinity.  Therefore,  a  fluctuation  in  acoustic  travel  time  is 
indicative  of  changes  in  the  environment  through  which  the 
acoustic  energy  has  passed.  Once  the  travel  time  fluctuations 
are  known,  inverse  mathematical  techniques  can  be  used  to 
infer  various  properties  of  the  ocean  medium.  [Ref.  1] 

One  source  of  travel  time  measurements  is  through  the  use 
of  explosive  devices.  Although  such  signals  are  easy  to 
generate,  they  may  not  provide  the  required  precision  owing  to 
dispersion  of  the  various  frequency  components  in  the 
impulsive  signal.  Additionally,  such  signals  are  difficult  or 
impossible  to  replicate.  [Ref.  2] 

A  better  method  that  has  been  employed  in  recent  years  is 
the  use  of  maximal -length  sequences .  The  technique  utilizes 
a  pseudorandom,  binary  sequence  as  the  phase  modulation  for  an 
electronically  generated  bandpass  signal.  Maximal-length 
sequences  are  well  suited  to  travel  time  measurements  by 


virtue   of   their   deterministic   nature,   autocorrelation 
properties  and  simplicity. 

The  goal  of  the  work  described  in  the  following  pages  is 
to  provide  sufficient  gain  to  permit  the  recovery  and  decoding 
necessary  for  accurate  travel  time  measurements  for  acoustic 
paths  of  extreme  length  (in  this  case  10,000  nautical  miles). 
Specifically,  the  programming  package  should  be  able  to: 


•  Exploit  the  spatial  structure  of  the  incident  wave  field 
to  discriminate  among  signals  arriving  from  different 
directions . 

•  Utilize  information  regarding  time  varying  array  tilt  and 
depth  to  estimate  the  array  geometry. 

•  Provide  a  stable  virtual  array  by  using  the  array  geometry 
and  modal  structure  of  the  immediate  environment. 

•  Distinguish  among  the  various  modal  components  of  the 
target  signal. 

•  Provide  both  time  domain  and  frequency  domain  analysis. 

•  Remain  sufficiently  flexible  so  as  to  permit  the 
processing  of  arrays  of  differing  construction  in  future 
tomography  experiments . 

•  Provide  some  measure  of  fault  tolerance  with  respect  to 
the  receiving  array. 


The  remainder  of  this  thesis  is  structured  as  follows: 

•  Chapter  II.  Acoustic  Wave  Propagation  Theory 

•  Chapter  III.  Modal  Beamforming 

•  Chapter  IV.  The  Heard  Island  Array 

•  Chapter  V.  Results  and  Conclusions 


Chapter  II  presents  an  introduction  to  acoustic  wave 
propagation.  The  approach  taken  is  that  of  normal  mode 
propagation  in  a  range  independent  channel. 

The  third  chapter  is  an  introduction  to  the  concept  of 
modal  beamforming.  It  reviews  frequency  domain  and  time 
domain  beamformers.  The  algorithm  utilized  in  this  study  is 
derived. 

Chapter  IV  describes  the  construction  and  subsequent 
deployment  of  the  receiving  array  used  in  this  experiment. 
Data  acquisition  and  preprocessing  of  the  acoustic  signal  are 
discussed. 

The  final  chapter  presents  a  description  of  the  array 
dynamics  encountered  and  frequency  domain  output  from  the 
software  developed.  Additionally,  proposals  for  system 
integration  and  performance  improvements  are  detailed. 

Appendix  A  contains  the  source  code  which  has  served  as 
the  body  of  this  work.  In  addition  to  the  beamformer,  various 
utility  programs  are  included  which  provide  for  array  geometry 
data  reduction. 

B.   THE  HEARD  ISLAND  EXPERIMENT 

The  performance  of  the  software  package  is  evaluated  on  a 
raw  data  set  (vice  synthetic  data).  The  data  was  acquired 
during  the  Heard  Island  Experiment  which  occurred  during  the 
period  January  23  -  February  2,  1991.  The  purpose  of  the 
experiment  was  to  determine  the  reliability  of  global  acoustic 


paths  for  tomographic  analysis.  Specifically,  it  is  desired 
to  establish  the  viability  of  these  transmission  paths  for  a 
proposed  multi-national  attempt  to  detect  a  decreasing  trend 
in  acoustic  travel  time.  Such  a  change  would  indicate  an 
overall  warming  along  the  path,  providing  the  first  convincing 
evidence  to  the  existence  of  global  warming  [Ref.  1]. 

The  signals  were  transmitted  from  the  vicinity  of  Heard 
Island  in  the  southern  Indian  Ocean.  The  site  was  selected 
because  it  is  central  to  a  web  of  open  water  paths  extending 
through  all  the  worlds  oceans  as  shown  in  Figure  1.1. 


Figure  1.1:  Heard  Island  Raypaths 


The  transmitting  ship  was  the  R/V  Cory  Chouest.  The 
transmitter  utilized  a  ten  element  vertical  array  with  a 
maximum  of  five  elements  active  at  any  time.  A  nominal  source 
level  of  209  dB  re  1  pPa  at  1  meter  was  realized.  Individual 
projectors  employed  were  HLF4LL  very  low  frequency  sources. 


A  variety  of  signal  types  were  sent,  including  continuous 
wave  and  maximal-length  sequences.  A  carrier  frequency  of  57 
Hz  was  chosen  both  for  its  low  absorption  losses  and  the 
ability  to  distinguish  it  from  the  50  and  60  Hz  frequencies 
generated  by  power  plants  world  wide. 

The  Monterey  component  of  the  Heard  Island  Experiment 
involved  a  collaboration  between  the  Naval  Postgraduate 
School,  the  Monterey  Bay  Aquarium  Research  Institute,  the 
Massachusetts  Institute  of  Technology,  Woods  Hole 
Oceanographic  Institution,  SAIC  and  the  Moss  Landing  Marine 
Laboratory.  The  R/V  Point  Sur  deployed  a  32  element  vertical 
array  approximately  70  nautical  miles  south-west  of  Monterey 
Bay. 


II.   ACOUSTIC  WAVE  PROPAGATION  THEORY 

A.   ACOUSTIC  PROPAGATION  IN  AN  OCEAN  WAVEGUIDE 
1.   The  Inhomogeneous  Wave  Equation 

The  homogeneous   linear  wave   equation    [Ref.  3], 

V2p=A^£,  (2.1) 

c2  at2 

incorporates  two  implicit  assumptions.  The  first  restricts 
the  sound  speed  to  a  constant  value  in  the  region  of  interest. 
Following  the  development  of  Coppens  [Ref.  4],  if  the  sound 
speed  is  now  permitted  to  vary  in  space  (but  not  time),  then 
c=c(x,y,z)  may  be  substituted  into  the  wave  equation  without 
loss  of  generality. 

The  second  implicit  assumption  restricts  the  use  of 
the  wave  equation  to  sound  fields  which  are  free  of  sources. 
Before  including  a  source  term,  the  selection  of  a  specific 
coordinate  system  is  appropriate. 

Given  that  sound  in  the  ocean  does  not  spread 
spherically  in  a  free  field,  but  radially  with  upper  and  lower 
boundaries,  implies  the  use  of  a  cylindrical  coordinate 
system.  To  further  simplify  the  problem  to  one  sufficient  for 
this  work,  the  sound  field  shall  be  assumed  to  exhibit  radial 
symmetry  about  the  source,  and  sound  speed  variations  shall  be 
restricted  to  depth,  c(z).    Based  on  these  restrictions, 


application  of  the  Lapacian  yields  the  homogeneous,  range 
independent  wave  equation, 


1£US\p+-^p  = 


i   a2 


(2.2) 


i  di\     drT      dz2  c2(z)    dt2 

The  inclusion  of  a  source  terra  requires  that  the 
inhomogeneous  wave  equation  reduce  to  homogeneous  form  in 
regions  which  are  free  of  sources  of  acoustic  energy.  This 
requirement  implies  that  the  delta  function  for  a  point  source 
located  at  r=0  and  z=zo  have  the  form  [Ref.  3] 


h{r-r, 


1 


6<r)6(z-z0), 

2  n  i  ° 


(2.3) 


where  r  is  the  spherical  radius  vector,  r  is  the  cylindrical 
radial  coordinate  and  z  is  the  cylindrical  depth  coordinate. 
Inclusion  of  this  term  in  the  wave  equation  yields  the 
Helmholtz  equation  for  a  monof requency  point  source  of  unity 
amplitude, 


±4-1x3- 

r  dr\     dz 


0) 


d2 
dz2     \c(z) 


p  =  --b(r)  b(z-z)  ejut 

r  ° 


(2.4) 
2 .   Solutions  to  the  Inhomogeneous  Wave  Equation 

Having  adopted  a  cylindrical  coordinate  system, 
equation  (2.4)  can  be  solved  by  separation  of  variables. 
Furthermore,  the  complete  solution  can  be  treated  as  a  linear 


combination  of  normal  modes, 

p(r,z,t)    =  e^^Rjr)  Zjz)  .  (2.5) 

m 

The  normal  modes  (ZJ  form  an  orthonormal  set  of 
eigenf unctions  in  z  which  satisfy  the  sourceless  Helmholtz 
equation, 


— +f-^-Kffl2>k  =  0,  (2.6) 

dz2     \c2(z)       m)    m 


subject  to  the  appropriate  surface  and  bottom  boundary 
conditions,  and  the  orthonormal  condition, 


fzn(z)Zjz)dz   =  8ffl(  (2.7) 


where  Km  is  the  eigenvalue  for  the  mth  eigenf  unction  (or  normal 
mode)  and  6„m  is  the  Kronecker  delta   function. 

'  nm 

A  closed  form  expression  for  the  Rm(r)  may  be  obtained 
by  substitution  of  equation  (2.6)  into  the  inhomogeneous  wave 
equation  (2.4),  multiplication  of  all  terms  by  Zn,  integrating 
over  z  and  application  of  the  normalization  condition  (2.7). 
The  result,  after  manipulation,  is 

\fr{r^y<K--*itU)Z^z0),  (2-8) 

the  inhomogeneous  Bessel's  equation.    It  has  the  known 
solution 

Rjr)    =  -jnZm(z0)H0^(Kmr)f  (2.9) 


where  H  <2)(k  r)  is  the  Hankel    function   of   the   second  kind  and 

o   v  m  ' 

order    zero.       Substitution  of  this  form  into  equation  (2.5) 
yields 

pd.Z.t)    =  -jne^^Zjz)  ZjzJ  H0<2)  (*„!)   .    (2.10) 

m 

Closed  form  solutions  for  the  Z  (z)  are  either 
difficult  or  impossible  for  all  but  select  cases. 
Fortunately,  efficient  numerical  algorithms  exist  which 
provide  for  the  rapid  solution  of  the  depth  dependent 
functions.  [Ref.  5] 


III.  MODAL  BEAMFORMING 

A.  PRELIMINARY  CONCEPTS 

Propagating  waves  can  be  modeled  as  functions  of  both 
space  and  time.  Conseguently,  the  attributes  of  such  models 
can  be  used  to  extract  information  from  real  wave  fields. 
Beamforming  exploits  the  temporal  and  spatial  characteristics 
of  a  specified  environment  to  enhance  a  particular  aspect  of 
the  wave  field  while  attenuating  undesirable  components.  Such 
an  operation  can  be  loosely  termed  constructive  reenf orcement 
of  the  desired  signal  or  noise  rejection.  The  information 
exploited  in  this  development  includes  the  modal  structure  of 
the  immediate  environment  and  the  instrument  data  supplied  by 
the  array  itself. 

B.  TIME  DOMAIN  BEAMFORMING 

The  acoustic  signal  received  at  a  given  hydrophone  can  be 
expressed  as 

p{t,x,y,z)    =  S(t,x,y,z)  +N{t,x,y,z)  ,  (3*1) 

where  p  is  the  pressure  at  the  hydrophone,  S  is  the  desired 
signal  and  N  is  the  local  noise  field.  All  are  functions  of 
time  and  hydrophone  location. 

If  one  now  considers  the  signal  received  at  an  array  of 
hydrophones,  eguation  (3.1)  takes  the  form 
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Pit)    =  Sit)  +N(t)  , 


(3.2) 


where  p,  S  and  N  are  all  N  x  1  vectors  representative  of  the 
individual  array  elements.  The  location  of  individual 
elements  is  implied  by  the  vector  index  for  a  given 
hydrophone . 

Classical  (or  plane  wave)  beamforming  incorporates  the  use 
of  a  complex  steering  vector  (E)  to  impose  a  phase  shift  on 
the  individual  elements  in  order  to  enhance  the  sensitivity  of 
the  array  to  signals  propagating  in  a  specific  direction, 


E  = 


J*N 


(3.3) 


The  output  signal  from  such  an  array  is 

bfpwit)    =  ET(S(t)  *mt))   •  (3.4) 

The  subscript  (pw)  indicates  that  the  beamformer  assumes  the 
presence  of  plane  waves.  Additionally,  the  implicit 
assumption  is  made  in  the  above  expression  that  the 
autocorrelation  function  of  the  noise  field  will  be 
identically  zero  for  any  two  (distinct)  hydrophones.  The 
success  or  failure  of  the  beamformer  will  be  a  direct 
reflection  of  the  validity  of  these  assumptions. 
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If  the  desired  signal  is  not  monof requency,  but  of  finite 
bandwidth,  simple  phase  shifting  may  not  be  suitable.  This  is 
particularly  true  if  the  target  signal  is  carrying  information 
such  that  the  communication  bit  time  is  comparable  to,  or 
greater  than  the  time  taken  for  the  wave  to  propagate  across 
the  array.  Under  such  circumstances,  the  beamformer  would 
impose  phase  shifts  in  distinct  communication  bit  segments 
simultaneously.  The  resulting  distortion  to  the  bandpass 
signal  is  unacceptable  in  most  communications  applications. 
Beamforming  under  these  conditions  requires  the  application  of 
true  time  delays  vice  phase  shifts. 

Assuming  an  array  of  spatial  extent  such  that  time  domain 
beamforming  is  required  implies  a  summation  of  the  form 

N 

bfpw(t)    =  £  (S(t-xn)  +N(t-xB))t  (3.5) 

n  =  l 

where  n  is  the  index  used  to  describe  the  array  elements  and 
rn  represents  the  steering  delay  applied  to  the  nth  element 
[Ref.  6].  This  arrangement  places  no  additional  restrictions 
on  the  target  signal  characteristics  or  array  construction  in 
a  given  physical  system. 

C.   THE  MODAL  BEAMFORMER 

The  signal  model  employed  above  is  incomplete,  in  that  it 
does  not  account  for  the  nonuniform  distribution  of  acoustic 
energy  in  the  sound  channel.    Vertical  arrays  inherently 
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sample  the  modal  structure  of  the  immediate  environment. 
Therefore,  one  needs  to  take  the  information  provided  by 
normal  mode  theory  into  account  when  the  acoustic  array  spans 
a  region  of  nonuniform  energy  density. 

The  appropriate  refinement  to  the  complex  steering  vector 
which  incorporates  the  modal  nature  of  the  signal  arrival  is 


E  = 


Zr{z2)ejif*   Z2(z2)eJ>" 


Zw(z1)eJ'*1" 
ZH(z2)ej*™ 


ZM(zN)ej*™ 


(3.6) 


Z.iz^e^"*  Z2(zN)ej*i* 

where  the  mth  column  of  the  matrix  is  the  generalized  steering 
vector  for  the  mh  normal  mode.  The  individual  steering 
vectors  now  include  amplitude  factors  (Zm(zn))  that  reflect  the 
value  of  the  eigenf unction  at  the  various  hydrophone  depths 
[Ref.  7].  The  nth  row  contains  the  values  of  those  steering 
vectors  at  the  depth  of  the  nth  hydrophone.  The  phase  terms 
(Tnm)  account  for  phase  shifts  required  for  the  beam  steering 
applicable  to  all  modes  (as  required  by  array  geometry)  and 
phase  delays  among  the  various  modes  due  to  their  individual 
propagation  speeds.   Specifically,  the  7nm  take  the  form 

ilr   =  6  +  a  ,  (3.7) 

where  0n  is  determined  by  the  hydrophone  location  and  am  is 
determined  by  the  differences  in  propagation  speeds  over  the 
transmission  path  for  the  normal  modes  present. 
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Modifying  the  steering  matrix  to  implement  time  delays 

(vice  phase  shifts)  is  straight  forward.   A  beamformer  could 

be  implemented  by  summing  the  output  of  individual  hydrophones 

(p(t))   at  this  point.    However,   to  do   so  implies  a 

deterministic  evaluation  of  the  a.   Such  an  evaluation  would 

n 

require  near  perfect  knowledge  of  the  conditions  over  the 
transmission  path  prior  to  commencement  of  the  experiment. 
Given  that  the  aim  of  tomography  is  the  mapping  of  these 
conditions,  such  an  approach  is  inappropriate. 

An  alternate  approach  is  to  utilize  the  columns  of  the 
matrix  individually  to  discriminate  among  the  modal  components 
of  the  signal.  Since  all  components  of  a  given  modal  steering 
vector  share  a  common  value  of  am,  these  phase  terms  may  be 
discarded  from  the  individual  elements  of  the  steering  vector. 
As  a  result  of  this  simplification,  the  modal  steering  vector 
is  completely  determined  by  the  Zm  (obtained  from  normal  mode 
modeling  of  the  immediate  environment)  and  the  desired  look 
direction  (determined  by  array  geometry) .  Having  established 
the  methodology,  and  replacing  the  phase  weights  with  steering 
delays,  the  output  of  the  modal   beamformer   is 

bfjt)    =  Y,Zm(zn)pnU-xn).  (3.8) 

n=l 

where  the  subscript  (m)  indicates  that  the  desired  output  is 
the  component  of  signal  representing  the  mth  normal  mode. 
Expanding   this   form   to   explicitly   include   the   terms 
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representing  the  signal  and  noise  components  of  the  acoustic 
field  yields 


N 


bfm(t)    =  £z-<za)Sa(£-Ta)+£z.<za)tfJI<t-tll)     (3.9) 


n=l  n-1 


where  Sn  and  Nn  represent  the  amplitudes  of  signal  and  noise 
at  the  nth  hydrophone.  The  delay  and  weighting  of  the  noise 
received  at  the  hydrophones  has  no  undesirable  effects.  If 
the  noise  between  any  two  hydrophones  was  uncorrelated  prior 
to  the  operation,  then  it  will  remain  so  after  weighting. 
Again,  the  degree  to  which  the  ambient  noise  field  is 
uncorrelated  will  ultimately  reflect  upon  the  performance  of 
the  beamformer. 
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IV.   THE  HEARD  ISLAND  VERTICAL  ARRAY 

A.   ARRAY  CONSTRUCTION 

The  array  deployed  for  this  experiment  consisted  of  32 
hydrophones  deployed  vertically,  each  having  a  nominal 
sensitivity  of  -170  dB  re  1  V/pPa.  The  element  spacing  was  45 
meters,  with  hydrophone  number  one  occupying  a  design  depth  of 
345  meters.  The  nominal  design  depth  for  hydrophone  number  32 
was  1740  meters. 

Instrumentation  on  the  array  consisted  of  two  sensors. 
The  upper  sensor  was  located  4.0  meters  above  hydrophone 
number  one.  It  recorded  ambient  pressure,  temperature,  tilt 
and  current  velocity.  The  lower  sensor  was  located  5.0  meters 
below  hydrophone  number  20.  It  recorded  tilt  (including 
direction),  pressure,  temperature  and  conductivity.  This 
sensor  appears  to  have  suffered  a  casualty  during  array 
deployment  which  rendered  its  tilt  data  suspect. 

Flotation  was  provided  by  one  main  syntactic  float  and  28 
syntactic  football  floats.  The  design  called  for  the  main 
float  to  reside  at  a  nominal  depth  of  230  meters. 
Approximately  half  of  the  football  floats  were  submerged. 
This  arrangement  rendered  the  array  neutrally  buoyant,  thus 
providing  a  measure  of  isolation  from  surface  wave  effects. 
The  array  was  directly  tethered  to  the  R/V  Point  Sur  with 


16 


floatation  devices  on  the  surface  portion.  Each  float  was 
equipped  with  a  flasher,  a  flag  and  a  radar  reflector.  The 
float  nearest  the  ship  had  the  addition  of  a  radio  beacon. 
The  procedure  called  for  the  ship  to  remain  approximately  1200 
meters  from  the  position  immediately  above  the  array.  During 
the  experiment  the  R/V  Point  Sur  was  configured  to  remain 
quiet  and  dead  in  the  water,  except  as  required  to  keep  clear 
of  the  array.  Figure  (3.1)  illustrates  the  general 
configuration . 


Main  float 


Upper  instrument 


Hydrophones 


Lower  instrument 


Ballast 


Figure  3.1:  Heard  Island  Receiving  Array 


17 


B.   INSTRUMENT  DATA 

Ideally  the  receiving  array,  when  deployed,  is  oriented 
vertically.  However,  due  to  currents,  the  requirement  to 
maintain  control  of  the  research  vessel  and  other  factors 
beyond  the  control  of  the  research  team,  the  array  is  often 
tilted  from  the  vertical.  Given  the  length  of  the  array,  a 
1.0°  tilt  results  in  a  24  meter  horizontal  displacement 
between  hydrophones  one  and  32.  This  displacement  is  nearly 
one  (carrier  frequency)  wavelength.  Additionally,  the 
direction  in  which  the  array  tilts  affects  the  component  of 
horizontal  displacement  which  is  collinear  to  the  propagation 
direction.  Horizontal  displacements  perpendicular  to  the 
direction  of  propagation  are  acceptable  if  one  assumes  that 
the  wave  field  is  fairly  uniform  is  this  direction. 
Therefore,  the  array  steering  utilized  in  this  work  seeks  only 
to  correct  for  the  component  of  displacement  in  the  direction 
of  signal  propagation. 

The  speed  with  which  the  wave  field  propagates  is  required 
prior  to  the  calculation  of  steering  delays.  The  appropriate 
speed  to  use  in  this  case  is  the  group  speed  of  the  carrier 
frequency  for  the  mode  in  question,  defined  as 


chc„ 


cgm=J^  (3.10) 


where  c   is  the  local  group  speed  for  the  mtn  normal  mode,  o> 
is  the  radian  frequency  and  Km  is  the  appropriate  eigenvalue. 
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The  amplitude  weight  applied  to  a  given  hydrophone  is  a 
function  of  the  mode  number  sought  and  hydrophone  depth.  The 
mode  shapes  (Zm(z))  are  required  prior  to  beamforming.  The 
numerical  algorithm  used  to  obtain  the  eigenvalues  and 
eigenf unctions  for  this  work  was  developed  by  Chiu  and  Ehret 
[Ref.  5].  It  utilizes  finite  differencing  of  the  propagation 
vectors  for  a  given  frequency  and  sound  speed  profile  to 
calculate  the  modal  structure  of  the  local  environment.  The 
sound  speed  profiles  were  measured  during  the  conduct  of  the 
experiment  by  the  research  team  aboard  the  R/V  Point  Sur. 
Figure  3.2  is  the  sound  speed  profile  used  for  this  analysis. 
Figures  3.3  through  3.5  are  the  output  of  the  normal  mode 
model.  Amplitude  weights  for  the  array  elements  were 
calculated  from  these  data  sets. 
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Sound  Speed  Profile 
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Figure  3.2:  Sound  speed  profile  at  receiving  array. 
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Figure  3.3:  Normal  modes  one  and  two. 
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Figure  3.4:  Normal  modes  three  and  four. 
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Mode  Five 
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Figure  3.5:  Normal  modes  five  and  six. 
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C.   DATA  ACQUISITION  AND  PREPROCESSING 

Once  deployed,  the  receiving  array  recorded  numerous  data 
sets,  each  lasting  in  excess  of  60  minutes.  Each  recording 
comprises  32  channels  of  acoustic  data.  The  hydrophone  output 
was  ported  to  a  patch  panel  and  preamplifier  where  a  fixed 
gain  of  26  dB  was  applied.  Following  this,  the  signal  passed 
through  a  variable  gain  amplifier  which  automatically  set  the 
signal  gain  to  optimize  the  10  volt  dynamic  range  of  the 
analog  to  digital  converter.  This  process  had  the  duel  effect 
of  enhancing  the  precision  of  the  data  and  preventing 
saturation  of  the  A/D  converters.  Subsequent  normalization  of 
the  data  set  to  hydrophone  output  levels  was  accomplished  by 
an  implementation  of  a  program  written  by  Keith  Von  der  Heydt 
of  Woods  Hole  Oceanographic  Institution. 

Analog  to  digital  conversion  utilized  a  sampling  frequency 
of  228  Hz.  The  signal  was  passed  through  an  analog  bandpass 
filter  prior  to  conversion.  This  prevented  aliasing  of 
frequencies  above  Nyquist  into  the  data  set.  Low  frequency 
noise  components  were  also  attenuated.  The  cutoff  frequencies 
for  the  bandpass  filter  were  10  and  80  Hz.  Data  storage 
utilized  both  optical  disk  and  magnetic  tape.  Data 
acquisition  and  preprocessing  of  the  data  sets  used  in  this 
work  was  conducted  by  Miller,  Chiu,  Frogner  and  others. 
Figure  3.6  illustrates  the  equipment  alignment.  [Ref.  8] 
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Figure  3.6:  Data  acquisition  system 
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V.   RESULTS  AND  CONCLUSIONS 

The  desired  result  of  this  research  is  the  development  of 
a  software  package  capable  of  beamforming  a  tomography  signal. 
The  beamformer  should  discriminate  among  the  various  modal 
components  of  the  signal,  thus  enabling  one  to  treat  the  modes 
individually.  Due  to  the  inherent  nature  of  vertical  arrays, 
the  software  needs  to  incorporate  information  pertaining  to 
the  time  varying  array  tilt  and  individual  hydrophone  depths. 
To  this  end,  the  beamformer  utilizes  a  duty  cycle.  Each  such 
cycle  commences  with  an  assessment  of  array  geometry.  This 
information  is  used  to  determine  the  appropriate  steering 
delays  and  hydrophone  weighting  coefficients.  Once 
established,  the  beamformer  processes  60  seconds  of  acoustic 
data  before  repeating  the  procedure.  In  this  manner,  the 
software  package  implements  a  quasi-stable  virtual  array  in 
which  the  virtual    aperture   remains  fixed  in  space. 

The  results  presented  examine  two  performance  aspects. 
The  first  pertains  to  array  geometry  description.  Calculation 
of  the  steering  delays  and  amplitude  weights  is  detailed. 
Additionally,  the  acoustic  output  of  the  array  is  compared  to 
that  of  a  single  hydrophone  on  the  sound  channel  axis,  thus 
allowing  for  a  quantitative  evaluation  of  the  ultimate 
performance  of  the  beamformer. 


26 


A.   HYDROPHONE  DELAY  AND  WEIGHTING 

Information  regarding  array  geometry  was  provided  by  two 
instrument  packages.  The  upper  sensor  was  located  4.0  meters 
above  hydrophone  number  one,  the  lower  sensor  5.0  meters  below 
hydrophone  number  20.  Figure  5.1  illustrates  the  instrument 
configuration  and  coordinate  system  used  in  determining  the 
array's  spatial  orientation. 

NORTH 


Upper  Instrument : 

NS  current  speed 
EW  current  speed 
Temperature 

Pressure 

Tilt 

Lower  Instrument: 

Temperature 
Pressure 
Conductivity 
NS  tilt 
EWtilt 


Figure  5.1:   Array  Instrumentation 
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1 .   The  Upper  Sensor 

Information  provided  by  the  upper  instrument  includes 
the  component  of  current  velocity  in  both  the  north-south  and 
east-west  directions,  ambient  pressure,  temperature  and  tilt. 
The  direction  in  which  the  array  tilts  is  determined  by  the 
current  direction.  The  primary  assumptions  being  that  the 
current  acting  on  the  array  is  independent  of  depth  and  that 
the  direction  in  which  the  array  tilts  is  completely 
determined  by  the  current  direction. 

In  making  these  assumptions,  one  considers  the  primary 
source  of  current  measured  by  the  instrument  to  be  a  result  of 
ship  positioning  adjustments.  As  mentioned  previously,  the 
R/V  Point  Sur  was  required  to  periodically  move  the  ship  in 
order  to  remain  clear  of  the  array.  Such  movements  were 
intended  to  straighten  the  tether  between  the  array  and  the 
towing  rig,  thus  preventing  the  tether  from  becoming  fouled  in 
the  ship's  propellers.  These  adjustments  resulted  in  a 
tensioning  of  the  tether.  It  is  this  tensioning  that  is 
considered  to  have  resulted  in  array  movement  and  consequent 
water  flow  across  the  instrument  package.  Based  on  this 
assessment,  the  assumption  that  tilt  direction  was  determined 
by  the  flow  measured  at  the  upper  instrument  is  consistent 
with  the  conditions  of  the  experiment.  Higher  order  models 
could  be  derived  for  array  dynamics,  however  this  research  has 
utilized  the  first  order  model  of  a  linear  array. 
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The  next  variable  pertinent  to  array  geometry  is  the 
pressure  measured  by  the  upper  sensor.  From  this  it  is 
straight  forward  to  calculate  the  sensor's  depth  using  the 
following  relationship, 

P - 1  a tmospheze  .5  .v 

99 

where  P  is  the  pressure   (in  pascals)   measured  by  the 

instrument,  p  is  the  density  of  seawater  and  g  is  the  force  of 

gravity. 

Figures  5.2  and  5.3  present  the  data  received  from  the 

upper  instrument  package  during  two  of  the  acoustic  recording 

periods.    The  recording  beginning  00:05  (GMT)  January  27 

coincides  with  the  reception  of  a  m-seguence  of  length  2047. 

A  continuous  wave  signal  was  received  during  the  recording 

beginning  15:05  (GMT)  January  27.   They  are  representative  of 

the  most  and  least  favorable  conditions  encountered  during  the 

course  of  the  experiment. 
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Figure    5.2:    Array   geometry    for   00:05    (GMT)    January   27 
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Figure  5.3:   Array  geometry  15:05  (GMT)  January  27. 
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2.   The  Lower  Sensor 

The  lower  sensor  was  intended  to  provide  a  redundant 
source  of  array  positioning  data.  Specifically,  it  measured 
tilt  in  the  north-south  and  east-west  directions,  temperature, 
pressure  and  conductivity.  As  discussed  previously,  the  lower 
instrument  appears  to  have  suffered  a  casualty  which  rendered 
the  tilt  data  suspect. 

Figure  5.4  indicates  that  an  event  occurred  which 
caused  the  sensor  to  provide  data  which  is  not  consistent  with 
the  assumed  array  geometry.  One  scenario  explaining  the  event 
is  that  the  sensor  failed  shortly  after  entering  the  water. 
Another  possible  explanation  is  that  the  array  became 
entangled  in  itself  upon  descent. 

All  array  tilt  data  was  stored  for  subsequent  analysis 
ashore.  Therefore,  the  research  team  had  no  knowledge  of  the 
apparent  failure  during  the  conduct  of  the  experiment.  As 
there  was  no  reason  to  suspect  the  instrument,  no 
investigation  was  made  to  determine  the  source  of  the 
inconsistencies . 
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Figure  5.4:   Lower  sensor  output  upon  deployment. 
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3.   Steering  Delays 

Time  delays  must  be  applied  to  each  of  the  hydrophones 
in  order  to  correct  for  array  tilt.  However,  only  the 
component  of  horizontal  displacement  collinear  to  the  signal 
propagation  direction  is  used  to  calculate  the  time  delay. 
This  reguirement  mandates  apriori  knowledge  of  the  propagation 
direction.  Prior  work  has  determined  that  the  expected 
direction  from  which  the  signals  would  arrive  is  217.0°  true 
[Ref.  9]. 

In  addition  to  array  geometry  and  propagation 
direction,  one  requires  a  nominal  propagation  speed  to  convert 
the  horizontal  displacements  to  time  delays.  The  appropriate 
speed  to  use  is  the  mode  group  speed,  as  previously  discussed. 
Group  speeds  for  this  study  were  obtained  from  eigenvalues 
calculated  by  the  normal  mode  model  courtesy  of  Chiu  and  Ehret 
[Ref.  5].  The  mode  group  speed  of  the  carrier  frequency  was 
used  to  calculate  time  delays.  Figure  5.5  presents  the  mode 
group  speeds  as  a  function  of  frequency  for  the  lowest  five 
modes . 
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Figure  5.5:   Mode  group  speeds. 
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4.   Hydrophone  Amplitude  Weighting 

Application  of  hydrophone  weights  utilizes  the  array 
geometry  to  determine  the  depth  of  each  array  element.  The 
value  of  the  depth  dependent  function  (Zm)  is  then  assigned  to 
the  appropriate  hydrophone  output.  This  manner  of  weighting 
enhances  the  array  sensitivity  to  the  mth  normal  mode. 

An  additional  consideration  encountered  in  hydrophone 
weighting  is  fault  tolerance .  Vertical  acoustic  arrays 
operate  in  a  hostile  environment,  thus  individual  hydrophones 
may  fail.  The  amplitude  weights  applied  to  the  array  must 
cope  with  the  failure  of  individual  elements.  To  this  end, 
the  beamformer  applies  a  weight  of  zero  to  elements  which  have 
been  evaluated  as  unreliable. 

Log  books  kept  by  the  research  team  aboard  the  R/V 
Point  Sur  indicating  suspect  hydrophones  were  verified  by  post 
processing  statistical  analysis  of  the  individual  elements. 
Several  hydrophones  were  found  to  have  failed  during  the 
course  of  the  experiment.  Specifically  elements  18  through  28 
were  evaluated  as  unreliable.  As  a  result  of  this 
determination,  only  hydrophones  one  through  17  were  used  to 
beamform  the  target  signals. 

B.   ACOUSTIC  PERFORMANCE 

The  preliminary  data  analysis  seeks  to  guantify  the 
beamformer ' s  performance  with  respect  to  a  single  hydrophone 
located  on  the  sound  channel  axis.    The  power  spectrum 
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calculated  from  the  output  of  hydrophone  five  was  compared  to 
the  spectra  calculated  from  the  beamformed  output  of  channels 
one  through  17.  The  array  gain    is  defined  as 


AG  =   I01og10 


SNR 


SNRn 


(5.2) 


where  SNRA  is  the  signal  to  noise  ratio  of  the  beamformed 
output  and  SNRn  is  the  signal  to  noise  ratio  of  the  nth 
hydrophone  [Ref.  10].  The  SNR  for  each  case  was  calculated  by 
integration  of  the  power  spectra  over  a  6.0  Hz  bandwidth 
centered  on  the  continuous  wave  frequency. 

The  array  gain  was  calculated  for  numerous  power  spectra, 
each  representing  60  seconds  of  acoustic  data.  Additionally, 
the  gain  was  calculated  for  the  lowest  five  normal  modes 
present.  Values  for  array  gain  range  from  a  low  of  2.5  dB  to 
a  high  of  8.5  dB .   The  nominal  array  gain  is   6  dB . 

Reasons  for  the  range  in  values  include  temporal 
variability  in  the  modal  components  of  the  signal,  errors 
introduced  in  the  steering  delays  and  errors  in  amplitude 
weighting.  Time  delay  and  amplitude  weight  errors  are  the 
direct  result  of  assumptions  made  regarding  array  geometry. 
The  primary  source  of  geometric  error  is  considered  to  involve 
the  assumption  that  tilt  direction  is  completely  determined  by 
current  direction.  Array  tilt  bearing  may  not  be  aligned  well 
with  current  direction  during  periods  when  the  direction  of 
flow  across  the  sensor  is  varying.   Therefore,  during  periods 
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of  variable  current  flow,  marginal  performance  can  be  expected 
of  the  beamformer. 

Figures  5.6  through  5.11  present  representative  power 
spectra  of  single  channel  data  and  beamformed  data  for  modes 
one  through  five  over  the  entire  range  of  frequencies 
analyzed. 
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Figure  5.6:   Hydrophone  five  power  spectrum. 
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Figure    5.7:      Mode   one  power   spectrum. 
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Figure    5.8:      Mode   two  power   spectrum. 
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Mode  Three  Power  Spectrum 
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Figure    5.9:      Mode   three  power   spectrum. 
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Figure    5.10:      Mode    four  power   spectrum. 
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Figure    5.11:      Mode    five  power   spectrum. 
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As  shown  in  the  previous  figures,  the  noise  field  at  the 
beamformer  is  not  incoherent.  This  fact,  coupled  with 
imperfect  geometric  array  description,  detracts  from  the 
calculated  array  gain. 

One  feature  worth  noting  is  the  decreasing  strength  of  the 
undesirable  tonal  components  with  increasing  mode  number, 
specifically  those  in  the  40  to  50  Hz  range.  Presumably, 
these  components  were  radiated  by  the  R/V  Point  Sur.  One 
hypothesis  is  that  the  strength  of  these  signals  in  mode  one 
is  the  result  of  grating  lobes  in  the  array's  beam  pattern  for 
this  mode.  The  amplitude  weights  applied  to  the  hydrophones 
are  similar  to  those  for  a  plane  wave  beamformer,  in  that 
there  is  no  phase  shift  applied  along  the  array  in  the  form  of 
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coefficients  of  differing  sign.  Therefore,  spatial  aliasing 
is  possible  for  frequencies  with  wavelengths  less  than  twice 
the  hydrophone  spacing.  For  the  Heard  Island  array,  this 
includes  frequencies  ranging  from  17  Hz  to  the  bandpass  cutoff 
frequency  of  80  Hz. 

As  the  mode  number  sought  increases,  these  frequency 
components  diminish  in  strength.  This  indicates  that  the 
rejection  of  coherent  noise  of  local  origin  improves  as  the 
number  of  180°  phase  changes  applied  to  the  hydrophone  output 
increases . 

Also  shown  in  these  figures  is  the  presence  of  distant 
sources  of  coherent  noise.  These  other  tonal  components  are 
seen  to  vary  both  as  a  function  of  mode  number  and  time.  This 
observed  behavior  is  consistent  with  modal  propagation  in  a 
range  dependent  environment.  Merchant  traffic  in  shallow 
water  and  offshore  drilling  activity  are  two  possible 
candidates  for  sources  of  coherent  noise  which  becomes  coupled 
into  the  sound  channel. 

Figures  5.12  through  5.17  present  the  same  power  spectra 
in  the  immediate  vicinity  of  the  carrier  frequency  (57  Hz). 
Again,  the  beamformed  output  indicates  successful  attenuation 
of  undesirable  components  of  a  coherent  noise  field, 
presumably  of  local  origin. 
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Figure    5.12:      Hydrophone    five   power   spectrum. 
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Figure  5.13:   Mode  one  power  spectrum. 
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Figure    5.14:      Mode   two  power   spectrum. 
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Figure  5.15:   Mode  three  power  spectrum. 
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Figure   5.16:      Mode    four  power   spectrum. 
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Figure  5.17:   Mode  five  power  spectrum. 
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C.  CONCLUSIONS 

Modal  beamf orming  is  a  valuable  technique  for  use  in  long 
range  tomography  experiments.  The  signals  presented  in  this 
work  traversed  the  longest  transmission  path  for  which 
reception  was  attempted  during  the  course  of  the  Heard  Island 
Experiment.  Although  no  travel  time  estimates  have  been  made, 
the  successful  detection  of  the  signals  indicates  that  large 
scale  tomography  experiments  are  feasible. 

The  primary  concern  of  this  work  was  the  development  of 
a  system  capable  of  detecting  the  acoustic  signals  emanating 
from  the  vicinity  of  Heard  Island.  This  has  been 
accomplished.  Additionally,  the  data  describing  the  complex 
dynamics  of  the  near  vertical  acoustic  array  has  implications 
lor  the  design  and  construction  of  vertical  arrays  for  use  in 
future  tomography  experiments.  The  most  troublesome  aspect  of 
this  work  has  been  the  accurate  determination  of  array 
geometry,  specifically  tilt  direction.  Describing  tilt 
direction  in  terms  of  measured  current  at  one  point  on  the 
array  is  considered  inadequate. 

D.  RECOMMENDATIONS 

Improvements  to  this  implementation  include  physical 
aspects  of  the  array  itself  and  performance  of  the  numerical 
algorithm.  Large  vertical  arrays  are  particularly  sensitive 
to  errors  in  estimated  geometry.  In  this  experiment,  an  error 
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of  0.25°  in  tilt  (assuming  perfect  knowledge  of  tilt 
direction)  creates  a  phase  error  of  45  degrees  between  the 
most  distant  elements.  This  error  is  more  destructive  to  the 
higher  modes  owing  to  their  greater  spatial  extent.  The 
problem  is  aggravated  by  guestionable  tilt  direction  data. 

Construction  of  future  vertical  arrays  should  incorporate 
multiple  instrument  packages  which  directly  measure  tilt  and 
direction.  As  this  was  the  case  with  the  lower  instrument  on 
the  Heard  Island  array,  it  is  most  unfortunate  that  it  failed. 
Additionally,  an  investigation  into  the  utility  of  higher 
order  models  describing  array  shape  would  be  illuminating. 
Incorporating  a  quadratic  model  for  array  shape  would  not 
impose  a  significant  computational  load  on  the  beamformer  and 
may  improve  the  estimated  position  of  individual  hydrophones. 

The  current  implementation  of  this  beamformer  utilizes 
third  order  polynomial  interpolation  during  the  application  of 
time  delays  to  individual  hydrophone  outputs.  As  such  it 
represents  the  greatest  computational  load  in  the  program. 
The  interpolator  is  based  on  Neville' s  algorithm  [Ref.  11]. 
No  attempt  was  made  to  minimize  execution  time.  As  a  result, 
the  software  requires  approximately  four  hours  to  process  one 
hour  of  acoustic  data.  A  more  efficient  interpolator  may 
reduce  execution  time  sufficiently  to  permit  this  technique  to 
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be  used  for  real  time  modal  beamforming.  To  this  end,  the 
software  is  assembled  from  modular  components,  thus  allowing 
for  improvements  on  a  function  by  function  basis. 
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APPENDIX  A 

The  following  programs  were  developed  during  the  course  of 
this  research.  They  include  the  modal  beamformer  and 
associated  utility  programs.  These  programs  may  be  obtained 
by  contacting  James  H.  Miller  at  his  address  in  the  initial 
distribution  statement. 

A.   THE  MODAL  BEAMFORMER 

1.   Operational  Considerations 

The  modal  beamformer  was  designed  with  portability  a 
prime  consideration.  As  a  result,  there  are  numerous  program 
parameters  which  must  be  set  by  the  end  user.  Among  these  are 
the  characteristics  of  the  physical  array  and  data  acquisition 
system.  All  such  user  selectable  parameters  are  contained 
solely  within  the  program  definitions.  The  parameters  appear 
below  in  the  same  order  in  which  they  appear  in  the  program. 


•  UNIX  VERSION:  selectively  compiles  either  unix  or  ansi 
compatible  code  blocks. 

•  ASCII:  selects  output  in  either  ascii  or  binary  format. 
Ascii  is  useful  for  data  export,  binary  saves  storage 
space . 

•  SIGNAL:  permits  user  to  enable  or  disable  the  beamformer' s 
time  series  output. 

•  SPECTRUM:  permits  the  user  to  enable  or  disable  the 
beamformer ' s  power  spectrum  output. 


49 


•  LOWER_SENSOR:  directs  the  program  to  load  or  ignore  lower 
tilt  sensor  data.  This  option  should  remain  OFF  as  there 
is  no  reliable  lower  sensor  data. 

•  VALIDATE:  enables  program  validation.  If  selected  the 
program  will  create  output  files  which  list  the  array 
geometry,  steering  delays,  hydrophone  weights  and 
calculated  mode  group  speed.  This  feature  is  useful  for 
program  debugging. 

•  ERROR_ESTIMATE:  if  selected,  this  feature  causes  the 
polynomial  interpolator  to  store  the  upper  bound  on  the 
mean  sguared  error  in  the  steering  delays.  This  feature 
doubles  the  output  of  the  beamformer. 

•  ON:  logical  switch  for  program  control. 

•  OFF:  logical  switch  for  program  control. 

•  INTERPOLATE:  selects  interpolator.  This  feature  permits 
the  installation  of  an  improved  interpolator  without  the 
requirement  to  search  for  and  replace  each  call  to  the 
function . 

•  ORDER:  indicates  the  degree  of  the  polynomial  used  in  the 
polynomial  interpolator. 

•  STEP:  indicates  the  number  of  points  on  either  side  of  a 
sequence  to  be  taken  when  estimating  derivatives. 

•  TINY:  prevents  division  by  zero  in  floating  point 
operations . 

•  PI:  used  for  trigonometric  recursions. 

•  RADIAN:  used  for  degree  to  radian  conversions. 

•  OFFSET:  indicates  the  distance  (in  meters)  between  the 
number  one  hydrophone  and  the  upper  tilt  sensor. 

•  DELTA_R:  indicates  array  element  spacing  (in  meters). 

•  CTD_OFFSET:  indicates  the  difference  in  depth  between  the 
sound  speed  profile's  first  data  point  and  the  depth 
increment  used  in  the  profile  (in  meters). 

•  SSP_LENGTH:  indicates  the  maximum  number  of  data  points 
that  an  eigenf unction  may  contain.  The  only  restriction 
applicable  to  this  parameter  is  available  memory. 
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•  EIGVAL_LENGTH:  indicates  the  maximum  number  of 
eigenvalues  that  the  program  will  accept.  The  only 
restriction  applicable  to  this  parameter  is  available 
memory. 

•  LOOK_DIRECTION:  the  compass  direction  (in  degrees  true) 
from  which  the  signal  arrives. 

•  TILT_BUFFER:  the  buffer  length  allocated  to  tilt  data. 
The  only  restriction  applicable  to  this  parameter  is 
available  memory. 

•  BUFFER_TIME:  the  length  of  time  (in  seconds)  represented 
by  the  acoustic  data  input  buffer.  This  period  represents 
the  output  buffer  length  plus  one  second  on  either  end  to 
permit  steering  delays  and/or  advances .  The  array  is 
currently  capable  of  end  firing  a  1500  meter  acoustic 
array.  This  parameter  must  be  an  integer.  Additionally, 
60  must  be  an  integer  multiple  of  (BUFFER_TIME-2 ) .  The 
array  has  a  duty  cycle  of  one  minute  between  design 
changes.  There  must  be  an  integer  number  of  output 
buffers  per  duty  cycle. 

•  F_SAMPLE :  sampling  rate  of  the  data  acguisition  system. 
This  value  must  be  an  integer. 

•  CHANNELS:   the  number  of  hydrophones  on  the  array. 

•  F_CARRIER:  the  carrier  frequency  of  the  target  signal. 
The  group  speed  used  in  calculating  steering  delays  is 
based  on  this  frequency.  This  must  be  entered  in  floating 
point . 

•  FFT_LENGTH:  the  length  of  the  FFT  used  if  frequency 
domain  output  is  selected.  This  value  must  be  a  power  of 
two  which  is  less  than  or  equal  to  the  number  of  samples 
in  a  single  channel  of  input  data. 

•  SWAP:  a  macro  used  in  calculating  a  FFT. 


The  following  program  functions  are  adaptations  of 
those  found  in  Reference  11: 

•  polint:  performs  the  polynomial  interpolation. 

•  realft:  calculates  Fast  Fourier  Transforms. 

•  fourl:   calculates  Fast  Fourier  Transforms 
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•  vector:  allocates  memory  for  one  dimensional  arrays. 

•  matrix:   allocates  memory  for  two  dimensional  arrays. 

•  free_matrix:    deallocates  memory  from  two  dimensional 
arrays . 

•  ExitOnError:   provides  abnormal  program  termination. 
2 .   Beamf ormer  Source  Code 


PROGRAM:  BEAMFORMER  vsn  1.0 
WRITTEN  BY:  Steven  Crocker 
LAST  UPDATE:  October  9,  1991 

This  program  takes  input  from  various  data  files  and  the  user.   It 
outputs  a  data  file.   The  inputs  are  a  number  of  channels  of  digital 
acoustic  data,  and  information  regarding  the  physical  characteristics 
and  geometry  of  the  receiving  array.   Additionally,  environmental  data 
in  the  form  of  normal  mode  eigenfunctions  and  eigenvalues  at  the 
receiving  array  are  required  to  operate  this  beamformer.   The  output 
is  a  single  channel  of  acoustic  data. 


#defi 
#defi 
#defi 
#defi 
fldefi 
#defi 
#defi 
#defi 
#defi 
#defi 
#defi 
#defi 
#defi 
#defi 
#defi 
#defi 
fldefi 
#defi 
#defi 
#defi 
#defi 
fldefi 
#defi 
#defi 
#defi 
#defi 
#defi 
#defi 


ne  UNIX  VERSION 

ne  ASCII  ON 

ne  BINARY  OFF 

ne  SIGNAL  ON 

ne  SPECTRUM  ON 

ne  L0WER_SENS0R  OFF 

ne  VALIDATE  OFF 

ne  ERROR_ESTIMATE  OFF 

ne  ON  1 

ne  OFF  0 

ne  INTERPOLATE  polint 

ne  ORDER  3 

ne  STEP  1 

ne  TINY  1.0e-25 

ne  PI  3.14159265359 

ne  RADIAN  57.2957795131 

ne  OFFSET  4.0 

ne  DELTA_R  45.0 

ne  CTD_0FFSET  0.5 

ne  SSP_LENGTH  2500 

ne  EIGVAL_LENGTH  230 

ne  L0OK_DIRECTI0N  217.0 

ne  TILT_BUFFER  120 

ne  BUFFER_TIME  12 

ne  F_SAMPLE  228 

ne  CHANNELS  32 

ne  F_CARRIER  57.0 

ne  FFT  LENGTH  2048 


/*  either  ANSI  or  UNIX  */ 

/*  select  output  mode  ON  or  OFF  */ 

/*  select  output  mode  ON  or  OFF  */ 

/*  either  ON  or  OFF  */ 

/*  either  ON  or  OFF  */ 

/*  either  ON  or  OFF  */ 

/*  either  ON  or  OFF  */ 

/*  either  ON  or  OFF  */ 

/*  logical  "switch"  */ 

/*  logical  "switch"  */ 

/*  polint  */ 
/*  order  of  interpolator         (odd)*/ 

/*  number  of  steps  for  derivatives  */ 

/*  prevents  division  by  zero     .  */ 

/*  for  freq  to  omega  conversions  */ 

/*  for  degree  to  radian  conversions  */ 

/*  dist  btwn  upper  sensor  and  phone  #1  */ 

/*  array  element  spacing  */ 

/*  diff  btwn  ctd  depth  inc  &  1st  depth  */ 

/*  max  number  of  pts  in  eigunfunction  */ 

/*  max  number  of  eigenvalues  */ 

/*  direction  from  which  signal  arrives  */ 

/*  max  length  of  tilt  data  vectors  */ 
/*  input  buffer  length  in  seconds  (int)*/ 
/*  sampling  frequency           (int)*/ 

/*  number  of  channels  processed  */ 

/*  carrier  frequency  */ 

/*  radix  2  <=  (BUFFER  TIME-2)*F  SAMPLE  */ 


fldefine  SWAP(a,b)  tempr=(a); (a)=(b); (b)=tempr 

#include<stdio.h> 
#  i  nc lude<ma I loc . h> 
#include<math.h> 

#if  defined  (  ANSI  ) 
#include<f loat.h> 
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#include<stdlib.h> 

int  getlnput(void); 

int  putOutput(void); 

int  processTi  It (float  **x,  float  **y,  float  **z); 

int  processModes( float  **z,  float  *weight,  float  *ptrC); 

int  dydx(float  *x,  float  *y,  float  *ddx,  int  points); 

int  polint(float  *xa,  float  *ya,  int  n,  float  x,  float  *y,  float  *dy); 

int  realft(float  *data,  int  n,  int  isign); 

int  fourKfloat  *data,  int  nn,  int  isign); 

float  *vector(int  length); 

float  **matrix(int  row,  int  col); 

int  free_matrix( float  **m,  int  row); 

void  ExitOnErroKchar  error_txt[]); 

tfelif  defined  (  UNIX  ) 

int   vectorO, 

■atrixO; 

getlnputO,  • 

putOutputO, 

processTi lt(), 

processModesC ), 

dydxO, 

polintO, 
realftO, 
fourK), 

f  ree_matrix(), 

ExitOnErrorO; 
#endif 

/*  Global  Variable  Declarations  */ 
float  **inSOUND,  *outSOUND,  *MeanSqError,  Max=0.0,  Min=1.0e25; 
int  lastMinute,  Minute=1,  f irstBuffer=1; 

mainO 

FILE  *fpVl; 

int  i,  j,  k,  n,  *shift; 

float  **x,  **y,  **z,  *indx,  *samples,  arg,  ans,  err,  *delay, 
*delta,  *weight,  *pwrSpectrum/  Cgroup; 


/*  Memory  Allocation  and  Tilt  Data  Processing  */ 
x=(f loat**)matrix( CHANNELS,  TILT_BUFFER); 
y=(float**)matrix(CHANNELS,  TILT_BUFFER); 
z=(f loat**)matrix(CHANNELS,  TILT_BUFFER); 
processTi  lt(x,  y,  z); 

if  (VALIDATE==ON) 

printf ("\nDumping  array  geometry  to  array. dat\n"); 

f pVl=fopen( "array. dat","wt" ); 

if (fpVl==NULL)  ExitOnErrorC'Error  opening  validation  file"); 

f printf (fpV1,"Channel\t\t  X\t\t  Y\t\t  Z\n"); 

for  (i=1;i<=lastMinute;i++) 
{ 

fprintfCfpVl,    "MINUTE:   %i\n",    i); 
for   (j=1;j<=CHANNELS;j++) 

fprintf(fpV1,"%i\t  %f\t  %f\t  5£f\n"/j/xCj]Ci3/yCj]Ci],zCj]Ci3); 
>  /*  for  */ 
fclose(fpVl); 
>  /*  if  */ 

/*  Memory  Allocation*/ 

weight= (float*) vector (CHANNELS); : 
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indx=(f loat*) vector (ORDER+1 ); 

delay= (float*) vector (CHANNELS); 

delta=(f loat*)vector(CHANNELS); 

outSOUND=(f  Loat*)vector( (BUFFER_TIME-1 )*F_SAMPLE); 

inSOUND=(float**)matrix(CHANNELS,  BUFFER_TIME*F_SAMPLE); 

if  ((shift=(int*)malloc((CHANNELS+1)*sizeof(int)))==NULL) 

ExitOnErrorC'Memory  allocation  failure  for  shiftC]."); 
for  (i=1;  i<=0RDER+1;  i++)  indx[i]=(f loat)i; 
if  (ERROR_ESTIMATE==ON) 

MeanSqError=(float*)vector((BUFFER_TIME-1)*F_SAMPLE); 
if  (SPECTRUM==ON)  pwrSpectrum=(f loat*)vector(FFT_LENGTH/2); 
whi  le(Minute<=lastMinute) 
{ 

/*  Mode  Data  Processing  */ 
processModesCz,  weight,  &Cgroup); 

/*  Calculate  delays,  shifts,  etc  */ 

for  (i=1;  i<=CHANNELS;  i++) 

{ 

delay[i]=x[i][Minute]/Cgroup;  /*  Time  delay   */ 

shift[i]=(int)(delay[i]*(float)F_SAMPLE);  /*  #   of  samples  */ 

/*  fraction  of  1  sample  */ 
deltaCi]=delayti]*(float)F_SAMPLE-(float)shiftCi]; 

>  /*  for  */ 

if  (VALIDATE==ON) 
{ 

printf ("Dumping  element  delays  to  delay. dat\n"); 

if  (f irstBuffer) 
{ 

if  ((fpVl=fopen("delay.dat",  "wt"))  ==  NULL) 

ExitOnErrorC'Error  opening  validation  file."); 
>  /*  if  */ 
else 
{ 

if  ((fpVl=fopen("delay.dat",  "at"))  ==  NULL) 

ExitOnErrorC'Error  opening  validation  file."); 

>  /•  else  */ 

fprintf(fpVl, "MINUTE:   %i\n",   Minute); 

fprintf (fpVl, 

"Channel\t  delay\t\t  int  shift\t  fraction  of  1  shift\n"); 

for  (i=1;i<=CHANNELS;i++) 
{ 

fpnntf(fpVl,   "%i\t  %e\t  %i\t%  e\n", 
i,delay[i3,shift[i],delta[i:); 

>  /*  for  */ 

fclose(fpVI); 

printf ("Dumping  phone  weights  and  grp  speed  to  modal. dat\n"); 

if(firstBuffer) 
{ 

if  ((fpVl=fopen("modal.dat",  "wt"))  ==  NULL) 

ExitOnErrorC'Error  opening  validation  file.Vn"); 

fprintf (fpVl, "Group  speed  for  F_CARRIER  is:  %g\n\n\n",Cgroup); 

>  /*  if  */ 
else 

{ 

if  ((fpV1=fopen("modal.dat","at"))  ==  NULL) 

ExitOnErrorC'Error  opening  validation  file.W); 
fprintf (fpV1, "\nHydrophone  weights  for  minute  %i\n", Minute); 

>  /*  else  if  */ 
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fprintf(fpV1,  "Channel  \tWeight\n"); 

for  (i=1;i<=CHANNELS;i++) 

fprintf(fpV1,  "z'.iVty.eNn",  i/UeightCiD^- 
fclose^pVl); 
>  /*  if  */ 

if  (SPECTRUM==ON) 

for  (k=1;k<=FFT_LENGTH;k++)  pwrSpectrum[k>0.0; 

for  (n=1;  n<=6G7(BUFFER_TIME-2);  n++) 
{ 

get  Input (); 

if (f irstBuffer)   /*  Produce  first  output  buffer  */ 
< 

for  (1=1;  i<=F_SAMPLE;  i++) 
{ 

outSOUND[i]=0.0; 

if  (ERROR  ESTIMATE==ON)  MeanSqErrorCi ]=0.0; 

>  /*  for  */  ~ 

for  (i=F  SAMPLE+1;  i<=(BUFFER_TIME-1)*F_SAMPLE;  i++) 
i 

outSOUND[i]=0.0; 

if  (ERROR_ESTIMATE==ON)  MeanSqError[i]=0.0; 

for  (j=1;  j<=CHANNELS;  j++) 
{ 

if  (delay[j]>=0.0) 
{ 

arg=(float)(0RDER+1)/2.0+(1-delta[j]); 
samples  =  &inS0UND[jHi-shift[j]-(0RDER+1 )/23; 
}  /*   if  */ 

else  if   (delayEjKO.O) 
{ 

arg=  (float  XORDER+D/2. 0+delta[j]; 
samples  =  &inS0UND[j][i-shift[j]-(0RDER+1)/2+1]; 
}  /*  else  if  */ 

INTERPOLATE(indx,  samples,  ORDER+1,  arg,  Sans,  &err); 
outSOUND[i]=outSOUND[i]+ans; 
if  (ERROR_ESTIMATE==ON) 

MeanSqError[i]=MeanSqError[i]+err*err; 
>  /*  for  */ 

if (fabs(outSOUNDCi])>Max)  Max=fabs(outSOUND[i]); 
if(fabs(outSOUND[i])<Min)  Min=fabs(outSOUND[i]); 
if  (ERROR_ESTIMATE==ON) 

MeanSqError[i]=MeanSqError[i]/(f loat) CHANNELS; 

>  /*  for  */ 
>  /*  if  */ 

else     /*  Produce  subsequent  output  buffers  */ 
{ 

for  (i=F_SAMPLE+1;  i<=(BUFFER_TIME-1)*F_SAMPLE;  i++) 
i 

outSOUNDCi-F_SAMPLE:=0.0; 

if  (ERROR_ESTIMATE==ON)  MeanSqError[i-F_SAMPLE]=0.0; 

for  (j=1;  j<=CHANNELS;  j++) 
{ 

if  (delay[j]>=0.0) 
{ 

arg=(float)(0RDER+1)/2.0+(1-deltaCj]); 

samples  =  &inS0UND[j][i-shift[j]-(0RDER+1)/2]; 
>  /*  if  */ 

else  if  (delay[j]<0.0) 
{ 

arg=  (float  >  (ORDER+1  )'/2 . 0+deltaCj]; 

samples  =  &inS0UND[j][i-shift[j]-(0RDER+1)/2+1]; 
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>  /*  else  if  */ 

INTERPOLATE(indx,  samples,  ORDER+1,  arg,  Sans,  Serr); 
out SOUND[ i-F_SAMPLE]=outSOUND[i-F_SAMPLE]+ans; 
if  (ERROR_ESTIMATE==ON) 

MeanSqErrorCi-F_SAMPLE]=MeanSqErrorCi-F_SAMPLE]+err*err; 
>  /*  for  */ 

if(fabs(outSOUNDCi-F_SAMPLE])>Max) 
Max=fabs(outSOUNDCi-F_SAMPLE:); 

if(fabs(outSOUND[i-F_SAMPLE])<Min) 
Min=fabs(out SOUNDC i-F_SAMPLE]); 
if  (ERROR_ESTIMATE==ON) 

MeanSqError:i-F_SAMPLE]=MeanSqError[i-F_SAMPLED/ 
(float) CHANNELS; 
>  /*  for  */ 
>  /*  else  */ 

if  (SIGNAL==ON)  putOutputO; 

if  (SPECTRUM==ON) 
{ 

window(outSOUND,  FFT_LENGTH); 

realftCoutSOUND,  FFT_LENGTH/2,  1); 

for  (k=0;k<FFT_LENGTH;k  +=  2) 

{ 
pwrSpectrum:k/2+13=pwrSpectrum:k/2+1]+ 
outS0UNDCk]*outS0UNDCk:+outS0UNDCk+1]*outS0UNDCk+1]; 

>  /*  for  */ 

>  /*  if  */ 

f irstBuffer=0;   /*  set  firstBuffer  to  false  */ 
>  /*  for  */ 

if  (SPECTRUM==ON) 
{ 

for  (k=1;  k<=FFT_LENGTH;  k++) 

pwrSpectrum:k]=pwrSpectrum:k]*(f loat)(BUFFER_TIME-2)/60.0; 

dumpSpectrum(pwrSpectrum); 
>  /*  if  */ 

pnntf("\t%i  minutes  of  input  data  processed. \n",  Minute); 
Minute++;       /*  increment  minute  counter  */ 
}  /*  while  */ 

printf ("EXECUTION  COMPLETE:  End  of  tilt  data  encountered\n"); 
printf ("Maximum  magnitude  encountered  was:  %e\n",Max); 
printf ("Minumum  magnitude  encountered  was:  %e\n",Min); 
exit(O); 

I irkick'k'k'k'kicic'k'k'k'kiricXicXXicicic'k'k'k    END    fn91  n   irkiric'k'k'kic'k'k'k'kic'k'k'k'k'k'k'k'kicicic'kirk'k'k'k / 
/icicicirk 

*  FUNCTION:  getlnputO 
* 

*  This  function  handles  all  acoustic  input.   It  also  provides  one  of 

*  two  normal  process  terminations  available  in  the  program.   (The  other 

*  is  located  in  mainO.) 


* 
* 

* 

Arguments: 

none 

* 

ir 

Return  value: 

0 

* 

Functions  called: 

vectorO 

ExitOnEr 

* 
* 
* 

it 

Definitions  called: 

ANSI 

F  SAMPLE 

BUFFER_TIME 

UNIX 
CHANNELS 

* 
* 
* 

Global  variables  ca 

lied: 

inSOUNDnn 
firstBuffer 

Min 
Max 
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*  Significant  memory  allocation:   diskBuf fer[] 

* 

A  A  A  A  A  J 

#if  defined  (  ANSI  ) 

int  getlnput(void) 

flelif  defined  (  UNIX  ) 

get  Input () 

tfendif 

{ 

int  i,  j,  buffer,  items; 

float  *diskBuffer; 

char  fileNameC80]; 

static  FILE  *fpStatic; 

if  (firstBuffer) 

printf ("Enter  file  name  for  input  acoustic  data:  "); 

if((scanf("Xs",  f i leName))==EOF) 

ExitOnErrorC'Fatal  error  in  scanfO"); 

printf ("\n\n"); 

if  ((fpStatic=fopen(fileName/  "rb"))  ==  NULL) 

ExitOnErrorC'Error  opening  INPUT  ACOUSTIC  data  file."); 

>  /*  if  */ 

/*  Memory  Allocation  */ 

if (firstBuffer)  buffer=BUFFER_TIME*F_SAMPLE*CHANNELS; 
else  buffer=(BUFFER_TIME-2)*F_SAMPLE*CHANNELS; 
diskBuffer=(f I oat*) vector (buffer ); 

items=f read((char*)(diskBuffer+1),  sizeof (float),  buffer,  fpStatic); 
if  (items==buffer)  /*continue*/; 
else  if (ferror(fpStatic)  !=  0) 

ExitOnErrorC'Error  encountered  while  reading  input  acoustic  data"); 
else  if(feof (fpStatic)  !=  0) 
{  ___   

pr  i  ntf  (  "  \n \ t icic%it&i(1c1rtrticicici(icic%ici:i(1:ici(ici(ieicicicicicicicicicicicicicitici(ic1cicic1e'kic'k  \n "  )  ." 

printf ("\t    End  of  File  reached:  EXECUTION  COMPLETED"); 

printf ("\t\t%i  minutes  of  data  processed\n",Minute-1); 

printf ("\t\t%i  bytes  of  data  discarded\n",  items*sizeof (float)); 

printf ("\t\tMaximum  magnitude  encountered  was:  %e\n",Max); 

printf ("\t\tMinimum  magnitude  encountered  was:  %e\n",Min); 

printf ("\t    End  of  File  reached:  EXECUTION  COMPLETED"); 

pr  i  nt  f  ( "\t*****^^********^*******^^  ) ; 

f close(fpStatic); 

exit(0); 

>  /*  else  if  */ 

else  ExitOnErrorC'Unknown  error  handling  acoustic  input  file."); 

for  (i=1;  i<=BUFFER_TIME*F_SAMPLE;  i++) 
i 

for  (j=1;  j<=CHANNELS;  j++) 
i 

if  (firstBuffer) 
{ 

inSOUND[j][i]  =  diskBuf fer[CHANNELS*(i-1)+j]; 

>  /*  if  */ 
else 

{ 

if(i<=2*F_SAMPLE) 

inS0UND[j][i]=inS0UNDCj][i+(BUFFER_TIME-2)*F_SAMPLE:; 
else 

inS0UND[j][i]=diskBufferCCHANNELS*(i-2*F_SAMPLE-1)+j]; 

>  /*  else  */ 
>  /*  for  */ 

>  /*  for  */ 
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/*  Deallocate  Memory  */ 
free((char*)diskBuffer); 
return(  0  ); 

/*****  END  getlnput  *****/ 
/AAA  AA**  AAAAAAAAAA4rAAA  AAA  A  / 

/AAAAA 

*  FUNCTION:  putOutputO 
* 

*  This  function  handles  all  acoustic  output.  Additionally,  it  outputs 

*  the  estimated  mean  squared  error  from  the  interpolators  (if  enabled). 


none 

0 

ExitOnErrorO 

ANSI 
F_SAMPLE 

OUtSOUNDCD 
f irstBuffer 


*  Arguments: 

* 

*  Return  value: 

* 

*  Functions  called: 
* 

*  Definitions  called: 

* 
* 

*  Global  variables  called: 
* 

* 

*  Significant  memory  allocation: 
* 

*if  defined  (  ANSI  ) 

int  putOutput(void) 

#elif  defined  (  UNIX  ) 

putOutputO 

#endif 

i 

int  i,  cut; 

char  f i leName[12],  mode[2D; 

static  FILE  *fpOutSound,  *fpMSE; 


if (f irstBuffer) 
t 

if  (ASCII==ON) 

{ 

modeC0]='W ; 
mode[1]='t ' ; 

>  /*  if  */ 

else  if  (BINARY==ON) 
{ 

mode[0]='w'  ; 

mode[1]='b'; 

>  /*  else  if  */ 

cut=1; 

printf ("Enter  file  name  for  output  data:  "); 

if((scanf("%s",  f i leName))==EOF) 

ExitOnErrorC'Fatal  error  in  scanfO"); 

printf("\n\n"); 

if  ((fpOutSound=fopen(fileName,  mode))  ==  NULL) 
ExitOnErrorC'Error  opening  OUTPUT  data  file."); 

if  (ERROR_ESTIMATE==ON) 
i 

printf ("Opening  file  error. dat\n\n"); 

if  ((fpMSE=fopen("error.dat",  mode))  ==  NULL) 
ExitOnErrorC'Error  opening  error.dat"); 

>  /*  if  */ 


UNIX 
BUFFER_TIME 

MeanSqErrorfJ 
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>  /*  if  */ 
else  cut=2; 

if  (ERROR_ESTIMATE==ON) 

if  (ASCII==ON) 
i 

for  (i=1;  i<=(BUFFER_TIME-cut)*F_SAMPLE;  i++) 
fprintf(fpMSE/"7.e\n"/  MeanSqErrorCi]); 

>  /*  if  */ 

else  if  (BINARY==ON) 
{ 

if(fwrite((char*)(MeanSqError+1)/sizeof(float)/(BUFFER_TIME-cut)* 

F_SAMPLE/fpMSE)==(unsigned)(BUFFER_TIME-cut)*F_SAMPLE)  ; 
else  if(ferror(fpMSE)  !=  0) 

ExitOnErrorC'Error  encountered  writing  error  data"); 
else  ExitOnErrorC'Unknown  error  handling  error  file."); 

>  /*  else  if  */ 

>  /*  if  */ 

if  (ASCII==0N) 
{ 

for  (i=1;i<=(BUFFER_TIME-cut)*F_SAMPLE;  i++) 
fprintf(fpOutSound,  "%e\n",  outSOUNDCi]); 

>  /*  if  */ 

else  if  (BINARY==0N) 
{ 

if(fwrite((char*)(outS0UND+1),sizeof(float)/(BUFFER_TIME-cut)* 

F_SAMPLE,fpOutSound)==(unsigned)(BUFFER_TIME-cut)*F_SAMPLE)  ; 
else  if (ferror(fpOutSound)  !=  0) 

ExitOnErrorC'Error  encountered  writing  output  acoustic  data"); 
else  ExitOnErrorC'Unknown  error  handling  acoustic  output  file."); 

>  /*  else  if  */ 
return(  0  ); 

>  /***iHnt*******************/ 

/*****  END  putOutput  *****/ 

/irkick'kirkirkickic'kirkirk'kickirk'k'k'k  / 
/irirk'k'k 

*  FUNCTION:  dumpSpectrumO 

* 

*  This  function  handles  dumps  the  signal  power  spectrum  (if  selected). 
* 


* 
* 

Arguments: 

pwrSpectrum 

* 

if 

Return  value: 

0 

• 

it 

Functions  called: 

ExitOnErrorO 

* 

Definitions  called: 

ANSI 

UNIX 

* 

F_SAMPLE 

FFTJ.ENGTH 

* 

Global  variables  Ce 

illed: 

none 

* 

Significant  memory 

allocat 

ion: 

none 

* 

/ 

#if  defined  (  ANSI  ) 

int  dumpSpectrum  (  float  *pwrSpectrum  ) 

#elif  defined  (  UNIX  ) 

dumpSpectrum  (  pwrSpectrum  ) 

float  *pwrSpectrum; 

#endif 

i 

int  i; 

static  float  sequenced .0; 

char  fileName[12],  mode [2]; 

static  FILE  *fp; 
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if (firstBuffer) 
{ 

if  (ASCII==ON) 
i 

modeCO]='w' ; 
modeC1]=,t'; 

>  /*  if  */ 

else  if  (BINARY==ON) 
{ 

modeCO]='W; 

mode[1]='b'; 

>  /*  else  if  */ 

printf ("Enter  file  name  for  output  SPECTRUM  data:  "); 

if((scanf("%s",  f i leName))==EOF) 

ExitOnErrorC'Fatal  error  in  scanfO"); 
printf("\n\n"); 

if  ((fp=fopen(fileName/  mode))  ==  NULL) 

ExitOnErrorC'Error  opening  OUTPUT  SPECTRUM  data  file."); 

printf ("Select  output  format:  \n"); 
printfC   Enter  0  for  MATLAB  compatible  output. \n"); 
printfC   Enter  1  for  GRAFTOOL  compatible  output. \n"); 
if((scanf("%i",  &i))==EOF) 

ExitOnErrorC'Fatal  error  in  scanfO"); 
if((i  !=  0)  &&  (i  !=  1))  ExitOnErrorC'Invalid  output  selection"); 

>  /*  if  */ 

if  (ASCII==0N) 
i 

if (firstBuffer  &S  format==1) 
i 

fprintf(fp,"0  "); 

for  (i=1,  i<=FFT  LENGTH/2;  i++) 

fprintf(fp,"%e  ", (float )(i-1 )*F_SAMPLE/FFT_LENGTH); 
fprintf(fp/"\n"); 

>  /*  if  */ 

if  (format==1)  fprintf (fp/'/ie  ",sequence); 

for  (i=1;i<=FFT_LENGTH/2;  i++) 

fprintf (fp,"%e  ",  pwrSpectrum[i]); 

fprintf (fp,"\n  "); 
sequence=sequence+1 .0; 
}  /*  if  */ 

else  if  (BINARY==0N) 
t 

if(fwrite((char*)(pwrSpectrum+1)/sizeof(float),FFT_LENGTH/2/fp)== 

(unsigned)FFT_LENGTH/2)  ; 
else  if (ferror(fpOutSound)  !=  0) 

ExitOnErrorC'Error  encountered  writing  output  SPECTRUM  data"); 
else  ExitOnErrorC'Unknown  error  handling  acoustic  SPECTRUM  file."); 

>  /*  else  if  */ 
return(  0  ); 

/*****  end  dumpSpectrum  ***•*/ 

*  FUNCTION:  processTi lt() 
* 

*  This  function  handles  all  array  tilt  data.   It  calculates  the  X,  Y,  Z 

*  coordinates  of  each  hydrophone  as  a  function  of  time.   The  coordinate 

*  system  is  oriented  such  that  X  points  toward  the  signal  "origin"  and 

*  Z  points  down. 
* 
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*  Arguments: 

* 


x[][] 

Z[][] 


y[][] 


Return  value: 
Functions  called: 

Definitions  called: 


ExitOnErrorO  matrixO 

vectorO  f  reejnatrixO 

ANSI  UNIX 

DELTA_R  LOWER_SENSOR 

RADIAN  CHANNELS 

LOOK_DIRECTION  OFFSET 
TILT  BUFFER 


lastMinute 


yy[][] 
tiltC] 
udepth[] 


*  Global  variables  called: 
* 

*  Significant  memory  allocation:   xx[][] 

*  zzllLl 

*  angled 

*  IdepthC] 
* 

A  A  A  AA / 

#if  defined  (  ANSI  ) 

int  processTi It (float  **x,  float  **y,  float  **z) 

#elif  defined  (  UNIX  ) 

processTi lt(x,y,z) 

float  **x,  **y,  **z; 

#endif 

< 

int  i,  j,  not  EOF; 

float  *tilt,  *angle,  *udepth,  *ldepth,  **Xx,  **yy,  **zz,  theta; 

char  f i  leName[12]; 

FILE  *fp1,  *fp2; 

/*  Open  Data  Fi  les  */ 

printf ("Enter  file  name  for  upper  tilt  sensor  data:  "); 

if((scanf("%s"/  f i  leName))==EOF) 

ExitOnErrorC'Fatal  error  in  scanfO"); 

printf("\n\n"); 

if  ((fp1=fopen(fileName,  "rt"))  ==  NULL) 

ExitOnErrorC'Error  opening  UPPER  TILT  data  file."); 

if  (LOWER_SENSOR==ON) 

printf ("Enter  file  name  for  lower  tilt  sensor  data:  "); 

if((scanf("%s",  f i leName))==EOF) 

ExitOnErrorC'Fatal  error  in  scanfO"); 

printf ("\n\n"); 

if  ((fp2=fopen(fileName/  "rt"))  ==  NULL) 

ExitOnErrorC'Error  opening  LOWER  TILT  data  file."); 
ldepth= (float*) vector (TILT_BUFFER); 
>  /*  if  */ 

/*  Memory  Allocation  */ 

tilt=(f loat*)vector(TILT_BUFFER); 

angle=(f loat*) vector (T1LT_BUFFER) ; 

udepth= (float*) vector (TILT_BUFFER); 

xx=(f loat** )matrix( CHANNELS,  TILT_BUFFER); 

yy=(f loat**)matrix( CHANNELS,  TILT_BUFFER); 

zz=(float**)matrix( CHANNELS,  TILT_BUFFER); 


1-1; 


/*  read  upper  tilt  data  */ 
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notE0F=1; 
while(notEOF) 
{ 
if(fscanf(fp1,M%g  %g  %g\n",&ti  ltr.i],&angle[i3/&udepth[i])    !=  EOF)   i++; 
else  notEOF=0; 
> 

if  (LOWER_SENSOR==ON) 
i 

j=1;  /*  read  Lower  tilt  data  */ 

notE0F=1; 

whiLe(notEOF) 

{ 

if(fscanf(fp2,  "%g\n",&ldepth[j])  !=  EOF)  j++; 
else  notEOF=0; 

>  /*  while  */ 

if  (i<=j)  lastMinute=i-1; 

else  lastMinute=j-1; 
> 
else  lastMinute=i-1; 

/************This  is  the  assumed  array  geometry:  LINEAR***********/ 

for  (j=1;  j<=CHANNELS;  j++) 

i 

for  (i=1;  i<=lastMinute;  i++) 
{. 

xx[j][i]=DELTA_R*(f loat)(j-1)*cos(angle[i]/RADIAN)* 

sin(tilt[i]/RADIAN); 
yy[j]CiD=DELTA_R*(f loat)( j-1 )*sin(angle[i]/RADIAN)* 

sin(tilt[i]/RADIAN); 
zz[j][i]=DELTA_R*(f loat)(j-1)*cos(tilt[i]/RADIAN)+ 
OFFSET* cos(tilt[i]/RADIAN)+udepth[i]; 

>  /*  for  */ 

>  /*  for  */ 

theta=(360.0-LOOK_DIRECTION)/RADIAN;  /*  coordinate  rotation  */ 

for  (j=1;  j<=CHANNELS;  j++) 

i 

for  (i=1;  i<=lastMinute;  i++)    /*  points  x  into  signal  */ 
i 

x[j][i]=xx[j][i]*cos(theta)-yy[j][i]*sin(theta); 
y[j]Ci]=xx[j][i]*sin(theta)+yy[j:[i]*cos(theta); 
zCj]Ci]=zz[j][i]; 

>  /*  for  */ 

>  /*  for  */ 

/*  Memory  Deallocation  */ 

if  (LOWER_SENS0R==0N)  f ree( (char*) Idepth); 

f ree((char*)tilt); 

f ree( (char*)angle); 

f ree( (char*)udepth); 

free_matrix(xx,  CHANNELS); 

free  matrix(yy,  CHANNELS); 

free_matrix(zz,  CHANNELS); 

fclose(fpl); 

if  (LOWER_SENSOR==ON)  fclose(fp2); 

returnC  0  ); 

/*****  END  processTilt  *****/ 

l~k  A  AAAAAAAAAAAAAAAAAAAAAAAfcA  / 
/*ArA"A"A 

*  FUNCTION:  processModesO 

* 

*  This  function  handles  the  normal  mode  data.   It  calculates  hydrophone 

*  weights  and  group  speed.   The  user  must  insure  that  the  depth  vector 

*  and  eigenfunction  vector  are  of  equal  length. 
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vectorO 
INTERPOLATEO 

ANSI 

PI 

EIGVAL_LENGTH 

CHANNELS 

Minute 


weightC] 


ExitOnErrorO 
dydxO 

UNIX 

SSP  LENGTH 

ORDER 

F  CARRIER 


Zn[] 

w[] 

dwdK 


*  Arguments:  z[][] 

*  ptrC 
* 

*  Return  value:  0 
* 

*  Functions  called: 

* 

*  Definitions  called: 
* 
* 
* 
* 

*  Global  variables  called: 
* 

*  Significant  memory  allocation:   depth[] 

*  OnOffC] 

*  Kr[] 
* 

ickickic  / 

#if  defined  (  ANSI  ) 

int  processModes(f loat  **z,  float  *weight,  float  *ptrC) 

#elif  defined  (  UNIX  ) 

processModes ( z , wei  ght , pt  r  C ) 

float  **z,  *weight,  *ptrC; 

flendif 

{ 

int  i,  j,  ptsEigVal,  set,  notEOF,  deadPhones,  weightNotAssigned; 

static  int  ptsEigFun; 

float  *w,  *Kr,  *dwdK,  err; 

static  float  depth[SSP_LENGTH+1],  Zn[SSP_LENGTH+1],  OnOf f [CHANNELS+1]; 

char  key,  f i  leName[12]; 

FILE  *fp1,  *fp2; 

if(firstBuffer==1) 
i 

w=(f loat*)vector(EIGVAL_LENGTH); 

Kr=(f loat*) vector (EIGVAL_LENGTH), • 

dwdK=(f loat*)vector(EIGVAL_LENGTH); 

printf ("Enter  file  name  for  normal  mode  data  (eigenfunction) :  "); 

i  f ( ( scanf ( "%s" ,  f  i leName) )==EOF ) 

ExitOnErrorC'Fatal  error  in  scanfO"); 

printf ("\n\n"); 

if  ((fp1=fopen(fi leName,  "rt"))  ==  NULL) 
ExitOnErrorC'Error  opening  NORMAL  MODE  data  file  (eigenfunction)"); 

printf ("Enter  file  name  for  normal  mode  data  (eigenvalues):  "); 

if ((scanf ("%s",  fi  leName) )==EOF) 

ExitOnErrorC'Fatal  error  in  scanfO"); 

printf ("\n\n"); 

if  ((fp2=fopen(fi leName,  "rt"))  ==  NULL) 
ExitOnErrorC'Error  opening  NORMAL  MODE  data  file  (eigenvalues)."); 

i=1;  /*  read  normal  mode  data  */ 

notE0F=1; 

while(notEOF) 

{ 

if(fscanf(fp1,  "%g  %g  \n",  &depth[i],  &Zn[i])  !=  EOF)  i++; 

else  notE0F=0; 
> 
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ptsEigFun=i-1; 

for(i=1;i<=ptsEigFun;i++)  depth[i]=depth[i]+CTD_OFFSET; 

1-1; 

notE0F=1; 

while(notEOF) 

{ 

if(fscanf(fp2,  "%g  %g  \n",  &w[i],  &Kr[i])  !=  EOF)  i++; 

else  notEOF=0; 
> 
ptsEigVal=i-1; 

for  (i=1;i<=CHANNELS;i++)  OnOff [i]=1 .0; 

printfC'Do  you  want  to  turn  off  any  hydrophones?  "); 

if((scanf("%s",&key))==EOF) 

ExitOnErrorC'Fatal  error  in  scanfO"); 

if(key=='y'  \ \    key=='Y') 
C 

printf ("\nHow  many  hydrophones  must  be  secured?  "); 

ifCCscanfC'/i",  &deadPhones))==EOF) 

ExitOnErrorC'Fatal  error  in  scanfO"); 

f or ( i  =1 ; i  <=deadPhones ; i  ++) 

printf ("\nEnter  hydrophone  number  to  secure:   "); 

if((scanf("%i",  8j))==E0F) 

ExitOnErrorC'Fatal  error  in  scanfO"); 

if  (j>CHANNELS  j ]  j<1) 

ExitOnErrorC'Bad  hydrophone  identification"); 
On0ff[j]=0.0; 
>  /*  for  */ 

>  /*  if  */ 

>  /*  if  */ 

for(i=1;i<=lastMinute;i++) 
{ 

if(depth[ptsEigFun]<zCCHANNELS][i]) 
{ 

printfC'Max  eigenfunction  depth  is:  %f \n",depth[ptsEigFun3); 

printfC'at  depth[]  index  of:  %i\n",ptsEigFun); 

printfC'Max  depth  of  phone  number  %i  is:  %f\n\n", 
CHANNELS, z[CHANNELS][i ]); 

ExitOnErrorC'Fatal  data  set  error"); 

>  /*  if  */ 

>  /*  for  */ 

for  (i=1;  i<=CHANNELS;  i++) 
i 

weightNotAssigned=1; 

while  ( j<=ptsEigFun  &S>  weightNotAssigned) 

£ 

if(zCi][Minute]<0.0  J|  depth[j]<0.0) 
i 

pnntf("i=%i\n",i); 

printf ("j=%i\n", j); 

printf ("Mi  nut e=%i\n", Minute) ; 

printf C'zCiDCMinute]   is:   %T\n",    z[i][MinuteD); 

printf C'depthCj]   is:   %f \n",depth[j]); 

printf ("Depth   less  than  zero  encountered  in  processModes. "); 

printf ("\n\n"); 
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ExitOnErrorC'Check  input  depths  for  coordinate  orientation"); 
>  /*  if  */ 


if(depth[j]<z[i][Mmute]  &&  depth[j+13>zCiKMinute]) 
{ 

set=(0RDER+1)/2; 

INTERPOLATE(&depth[j-setD,   &Zn[j-set],   ORDER+1,   zCiDCMinute], 
ftweightCi],   &err); 

weightNotAssignedO; 

>  /*  if  */ 

>  /*  while  */ 

>  /*  for  */ 

for  (i=1;  i<=CHANNELS;  i++)  weight[i]=On0ff [i]*weight[i]; 

if  (Minute==1) 

dydxCKr^dwdl^ptsEigVal); 
for  (i=1;i<=ptsEigVal;i++) 
{ 

if  (w[i]<2.0*PI*F_CARRIER  &&  w[i+1 ]>2.0*PI*F_CARRIER) 
{ 

set=(0RDER+1)/2; 

INTERPOLATE(&w[i -set], &dwdK[i-set], ORDER+1, 
2.0*PI*F_CARRIER,ptrC,&err); 

>  /*  if  */ 
}  /*  for  */ 

f ree((char*)w); 

free((char*)Kr); 

free((char*)dwdK); 

>  /*  if  */ 
returnC  0  ); 

y    /AAAAAAAAAAAAAAA^AAAAAAAAAAA'A'/ 

/*****  END  processModes  *****/ 

/irk  ttirirkirkickirkirkicirkicicirkirirkiricic  / 
/ickirtrk 

*  FUNCTION:  dydxO 
* 

*  This  function  estimates  derivatives. 
* 

*  Arguments: 
* 
* 

*  Return  value: 
* 

*  Functions  called: 
* 

*  Definitions  called: 
* 
• 

*  Global  variables  called:        none 
* 

*  Significant  memory  allocation:   none 
* 

#if  defined  (  ANSI  ) 

int  dydx(float  *x,  float  *y,  float  *ddx,  int  points) 

#elif  defined  (  UNIX  ) 

dydx(x,y,ddx, points) 

float  *x,  *y,  *ddx; 

int  points; 

#endif 

{ 

int  n; 


xc: 

y" 

ddx[] 

points 

0 

ExitOnErrorO 

ANSI 

UNIX 

STEP 
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for  (n=1;n<=points;n++) 
i 

if  ((n>=STEP)  &&  (n<=points-STEP))  /*center*/ 

ddx[nD=(y[n+STEP]-y[n-STEP])/(xCn+STEP]-xCn-STEP]); 

else  if  (n<STEP)  /*beginning*/ 

ddx[n]=(y[n+STEP3-y[1])/(x[n+STEP]-x[1]); 

else  if  (n>points-STEP)  /*end*/ 

ddx[n]=(y[points:-yCn-STEP])/(x[points:-x[n-STEP]); 

else 

ExitOnErrorC'Index  error  in  dydx");   /*  sanity  check  */ 
>  /*  for  */ 
return(  0  ); 

/*****  end  dydx  aaaaa/ 

/A  A  A  A  A 

*  FUNCTION:    polintO 
* 

*  This  function  performs  polynomial   interpolation. 


* 
* 
* 
* 

Arguments: 

xa[] 
n 

y 

ya[] 

X 

dy 

* 
* 
* 

Return  value: 

0 

Functions  called: 

vectorO 

ExitOnErrorO 

* 
* 
* 
* 
* 

Definitions  called: 

ANSI 

UNIX 

Global  variables  called: 

none 

Significant  memory  allocation: 

d[] 

c[] 

icirtckic  / 

tfif  defined  (  ANSI  ) 

int  polintC  float  *xa,  float  *ya,  int  n,  float  x,  float  *y,  float  *dy) 

#elif  defined  (  UNIX  ) 

polint(xa,ya,n,x,y,dy) 

float  *xa,  *ya,  x,  *y,  *dy; 

int  n; 

#endif 

{ 

int  i,  m,  ns=1; 

float  den,  dif,  dift,  ho,  hp,  w; 

float  *c,  *d; 

dif=fabs(x-xa[1]); 

c=(f loat*) vector (n), • 
d=(f loat*) vector (n); 

for  (i=1;  i<=n;  i++) 
{ 

if  ((dift=fabs(x-xa[i]))  <  dif) 

{ 

ns=i; 
dif=dift; 

>  /*  if  */ 

c[iD=yaCi]; 

d[i3=ya[i]; 
>  /*  for  */ 
*y=ya[ns— ]; 
for  (m=1;  m<n;  m++) 
{ 
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for  (1=1;  i<=n-m;  i++) 
{ 

ho=xa[iD-x; 
hp=xa[i+m]-x; 
w=c[i+1]-d[i]; 
if  Uden=ho-hp)==0.0) 

ExitOnErrorC'Error  in  routine  POLINT"); 
den=w/den; 
dCi]=hp*den; 
c[i]=ho*den; 
>  /*  for  */ 

*y  +=  (*dy=(2*ns  <  (n-m)  ?  cCns+1]  :  d[ns— ])); 
>  /*  for  */ 
f ree((char*)d); 
f ree((char*)c); 
return(  0  ); 

/*****  END  polint  ****/ 

/III     LJilA    1     I     I     I     I   liiliLii  JiJnl.JiX     / 
Iflf  JtltXltTXXXrtTXxxXXTxxX  / 

/irtckick 

*  FUNCTION:  windowO 
* 

*  This  function  applies  a  Blackman  window  to  a  vector. 
* 


*  Arguments: 

* 

Return  value: 


* 

* 

* 

* 

* 

* 

* 

* 

* 

*  Significant  memory  allocation: 

* 

"kickicic  / 


Functions  called: 
Definitions  called: 

Global  variables  called: 


dataCD 

0 

none 

ANSI 
PI 

none 

none 


N 


UNIX 


#if  defined  (  ANSI  ) 

int  window(float  *data,  int  N); 

flelif  defined  (  UNIX  ) 

windowC  data,  N  ) 

float  *data; 

int  N; 

tfendif 

{ 

int  n; 

for  (n=0;  n<N;  n++) 
i 

data[n+1]=data[n+1D*(0.42+0.5*cos(2.0*PI*(float)(n-N/2)/(float)(N-1)) 
+0.08*cos(4.0*PI*(f loat)(n-N/2) /(float) (N-1))); 
>  /*  for  */ 
return  (  0  ); 


/*****  end  window  *****/ 


*  FUNCTION:  realftO 
* 

*  This  function  calculates  FFT's 
* 

*  Arguments: 
* 

* 

*  Return  value: 
* 


dataC: 
i  s  i  gn 
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*  Functions  called:  fourK) 
• 

*  Definitions  called:  ANSI  UNIX 

*  Global  variables  called:        none 
* 

*  Significant  memory  allocation:   none 

* 

#if  defined  (  ANSI  ) 

int  realft(float  *data,  int  n,  int  isign) 

tfelif  defined  (  UNIX  ) 

realft(data,  n,  isign) 

float  *data; 

int  n,  isign; 

flendif 

<. 

int  i,  i1,  i2,  i3,  i4,  n2p3; 

float  d=0.5,  c2,  h1r,  Mi,  h2r,  h2i; 

double  wr,  wi,  wpr,  wpi,  wtemp,  theta; 

theta=3.141592653589793/(double)n; 

if  (isign==1) 

{ 

c2  =  -0.5; 

fourKdata,  n,  1); 

>  /*  if  */ 
else 

i 

c2=0.5; 

theta  =  -theta; 

>  /*  else  */ 
wtemp=sin(0.5*theta); 
wpr  =  -2.0*wtemp*wtemp; 
wpi=sin(theta); 

wr=1 .0+wpr; 

wi=wpi; 

n2p3=2*n+3; 

for  (i=2;  i<=n/2;  i++) 

{ 

i4=1+(i3=n2p3-(i2=1+(i1=i+i-1))); 

h1r=c1*(dataCi1]+dataCi33); 

h1i=c1*<dataCi2]-data[i^]); 

h2r  =  -c2*(data[i2]+data[iA3); 

h2i=c2*(data[i1]-data[i3]); 

dataC  i1 ]=h1 r+wr*h2r-wi*h2i ; 

dataCi2]=h1 i+wr*h2i+wi*h2r; 

data[i3]=h1r-wr*h2r+wi*h2i; 

dataEiA]  =  -h1i+wr*h2i+wi*h2r; 

wr=(wtemp=wr)*wpr-wi*wpi+wr; 

wi=wi*wpr+wtemp*wpi+wi; 

>  /*  for  */ 
if  (isign==1) 

<. 

datad ]=(h1 r=data[1 ] )+dataC2]; 
data[2]=h1r-data[2]; 

>  /*   if  */ 
else 

{ 

dataH  ]=d*(  (hi  r=data[1  ]  )+data[2]  ); 

data[2]=c1*(h1r-data[2]); 

four1(data,n,-1); 

>  /*  else  */ 
return  (  0  ); 

~y    /iririrkirkirkirirkirkickirkicirirk'k  / 

/*+•**  end  realft  *****/ 
/icicirk'k'frfrkirfrk'k'kiciric'X+'klrk'k  / 

/* 
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*  FUNCTION:   fourlO 

* 

*  This  function  calculates  FFT's 
* 

*  Arguments:  data[] 

*  isign 
* 

*  Return  value:  0 
* 

*  Functions  called:  none 
* 

*  Definitions  called:  ANSI 

*  SWAP 
* 

*  Global  variables  called:        none 
* 

*  Significant  memory  allocation:   none 
* 


nn 


UNIX 


/ 

#if  defined  (  ANSI  ) 

int  fourKfloat  *data,  int  nn,  int  isign) 
flelif  defined  (  UNIX  ) 
fourKdata,  nn,  isign) 
float  *data; 
int  nn,  isign; 
#endif 
< 

int  n,  mmax,  m,  j,  istep,  i; 

double  wtemp,  wr,  wpr,  wpi,  wi,  theta; 

float  tempr,  tempi; 

n=nn  «  1; 

for  (i=1;  i<n;  i+=2) 

<. 

if  (j  >  1) 
{ 

SWAP(data[j3,data[i]); 

SWAP(data[j+1],data[i+1]); 

>  /*  for  */ 
m=n  »  1; 

while  (m  >=  2  &&  j  >  m) 
{. 

j  -=  m; 

m  »=1; 

>  /*  while  */ 
3  +=  m; 

>  /*  for  */ 

mmax=2; 

whi le  (n  >  mmax) 

i 

istep=2*mmax; 

theta=6.28318530717959/(isign*mmax); 

wtemp=sin(0.5*theta); 

wpr  =  -2.0*wtemp*wtemp; 

wpi=sin(theta); 

wr=1 .0; 

wi=0.0; 

for  (m=1 ;  m<mmax;  m+=2) 

i 

for  (i=m;  i<=n;  i+=istep) 

j=i+mmax; 

tempr=wr*data[j]-wi*data[j+1]; 
tempi=wr*dataCj+1]+wi*data[j]; 
data[j]=dataCi!l-tempr; 
data[ j+1]=data[ i+1] -tempi ; 
data[i]  +=  tempr; 


69 


data[i+1]  +=  tempi; 
>  /*  for  */ 

wr=(wtemp=wr)*wpr-wi*wpi+wr; 
wi=wi*wpr+wtemp*wpi+wi; 
>  /*  for  */ 
mmax=istep; 
>  /*  while  */ 
return   (  0  ); 

/*****  end  fourl  *****/ 

/irtckick 

*  FUNCTION:  vectorO 
* 

*  This  function  allocates  memory  for  UNIT  OFFSET  vectors. 
* 

*  Arguments:  length 

*  Return  value:  *v 
* 

*  Functions  called:  ExitOnErrorO 
* 

*  Definitions  called:  ANSI  UNIX 
* 

*  Global  variables  called:        none 

*  Significant  memory  allocation:   v[] 
* 

irk~kirk  / 

#if  defined  (  ANSI  ) 

float  *vector(int  length) 

#elif  defined  (  UNIX  ) 

vector(  length) 

int  length; 

#endif 

{ 

float  *v; 

i f  ( (v=(f loat*)mal loc(  ( length+1 )*si  zeof (f loat ) ) )==NULL) 
ExitOnErrorC'Memory  allocation  failure  in  vectorO."); 
#if  defined  (  ANSI  ) 

return  v; 
#elif  defined  (  UNIX  ) 

return  (long  mt)v; 
flendif 

y    /irkickirkickirkickirkicirk'k'k'k'k'k  j 

/*•***  END  vector  *****/ 

/  A  A  A  A  A  A  kk"kk"ki(  A  AAA  kk"kk"kk  / 
/A  A  A  A  A 

*  FUNCTION:  matrixO 
* 

*  This  function  allocates  memory  for  UNIT  OFFSET  2-D  arrays. 
* 

*  Arguments:  row  col 
* 

*  Return  value:  **m 
* 

*  Functions  called:  ExitOnErrorO 
* 

*  Definitions  called:  ANSI  UNIX 
* 

*  Global  variables  called:        none 
* 

*  Significant  memory  allocation:   m[][] 
* 

A  A  A  A  A  / 

#if  defined  (  ANSI  ) 

float  **matrix(int  row,  int  col) 
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row 


flelif  defined  (  UNIX  ) 

matrix(row,col) 

int  row,  col; 

#endif 

{ 

int  i; 

float  **m; 

if  ((m=(float**)malloc( (unsigned) (row+1)*sizeof (float*) ))==NULL) 

ExitOnErrorC'Al location  failure  1  in  matrixO"); 
for  (i=1;  i<=row;  i++) 
i 

if  ((mCi]=(float*)malloc( (unsigned) (col+1)*sizeof (float )))==NULL) 
ExitOnErrorC'Allocation  failure  2  in  matrixO"); 
>  /*  for  */ 
#if  defined  (  ANSI  ) 

return  m; 
#elif  defined  (  UNIX  ) 
return  (long  int)m; 
flendif 

/*****  END  matrix  *****/ 
I irkirkirkirk^'k'k~kk"kirkirkk"kk"k  / 
/kkkkk 

*  FUNCTION:  f reejnatrixO 
• 

*  This  function  deallocates  memory  from  UNIT  OFFSET  2-D  arrays. 
* 

*  Arguments:  m[][] 
* 

*  Return  value:  0 
* 

*  Functions  called:  none 
* 

*  Definitions  called:  ANSI 
* 

*  Global  variables  called:        none 
* 

*  Significant  memory  allocation:   none 
* 

k  k  kirk/ 

#if  defined  (  ANSI  ) 
void  f ree_matrix(f loat  **m,  int  row) 
#elif  defined  (  UNIX  ) 
f ree_matrix(m,row) 
float  **m; 
int  row; 
#endif 
{. 

int  i; 

for(i=row;  i>=1;  i — ) 
f ree((char*)m[i]); 
free((char*)m); 
return(  0  ); 

y    /kkkkk  kk  kkkkkkk  k  k  kkirkk"  k  k  k  k  k  k  / 

/*****  END  free  matrix  *****/ 
/icicirkicirlrkiciclrk'kirklrk^'kicklciricickic  / 


UNIX 


/ 

*  FUNCTION:  ExitOnErrorO 
* 

*  This  function  performs  an  abnormal  process  termination. 
* 

*  Arguments:  error_txt[] 
* 

*  Return  value:  none 
* 

*  Functions  called:  none 
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*  Definitions  called:  ANSI 

*  Global  variables  called:        none 
* 

*  Significant  memory  allocation:   none 
* 

A  A  A  A  A  / 

#if  defined  (  ANSI  ) 

void  ExitOnErrorCchar  error_txt[]) 

#elif  defined  (  UNIX  ) 

ExitOnError(error_txt) 

char  error_txt[]; 

flendif 

{ 

fprintf (stderr, "Program  run-time  error  ...\n"); 

fprintf (stderr,"%s\n",error_txt); 

fprintf (stderr,". . .now  exiting  to  system. .. \n"); 

exit(O); 

/*-*■***  END  ExitOnError  AAAAA/ 

/^kick'k'kic'kXic'kkkkickirkk'kkk'k'k'kkk  / 


UNIX 


B .       DECOMMA . C 

This  program  removes  commas  and  colons  from  the  output  of 
the  upper  tilt  sensor  of  the  Heard  Island  Array. 


/irkirirk 

*     PROGRAM  DECOMMA. C  vsn  1.0 
WRITTEN  BY:  Steven  Crocker 
LAST  UPDATE:  August  3,  1991 


* 
* 
* 
* 
* 

A AAA  A  / 


This  program  takes  any  ASCII  file  and  copies  it  to  a  user 
defined  file.   It  removes  the  commas  and  colons. 


/*  either  ANSI  or  UNIX  */ 
/*  max  length  for  file  names  */ 
/*  ASCII  codes  */ 


fldefine  VERSION  ANSI 
fldefine  C  25 
#define  COMMA  44 
^define  COLON  58 
^define  SPACE  32 
#include<stdio.h> 
//include  <malloc.h> 


#if  (VERSION==ANSI) 

#include<stdlib.h> 

void  ExitOnErrorCchar  error_txt[]); 

tfendif 

mainO 
{ 

FILE  *sacm_fp,  *out_fp; 

char  c,  sacm_f i leCC],  out_f i leCC]; 

printf ("\n\n\nEnter  file  name  for  input  SACM  data.Xn"); 

scanf ("%s",sacm_f i le); 

printf ("XnXn"); 

printf ("Enter  file  name  for  output  SACM  data.Xn"); 

scanf ("%s",out_f i le); 

if  ((sacm_fp=fopen(sacm_f ile,  "rt"))  !=  NULL)  { 

if  ((out_fp=fopen(out_f ile,  "wt"))  !=  NULL)  -C 

whi  le((c=fgetc(  sacm_fp  ))  !=  EOF)  -C 
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if((c==COMMA)  |i  (c==COLON))  c=SPACE; 
fputc(c,  out_fp); 
>  /*  while  */ 
>  /*  if  */ 

else  ExitOnErrorC'Error  opening  output  data  file"); 
}  /*  if  */ 

else  ExitOnErrorC'Error  opening  input  data  file"); 
f close(sacm_fp); 
fclose(outjfp); 
exit(O); 


#if  (VERSION==ANSI) 

void  ExitOnErrorCchar  error_txt[]) 

#elif  (VERSION==UNIX) 

ExitOnError(error_txt) 

char  error_txtCD; 

#endif 

{ 

fprintf (stderr, "Program  run-time  error  ...\n"); 

fprintf (stderr,"%s\n"/error_txt); 

fprintf (stderr,". . .now  exiting  to  system. .. \n"); 

ex  i  t ( 1 ) ; 
> 


C .   SACM1 . C 

This  program  reads  the  data  output  by  the  upper  tilt 
sensor  and  formats  it  for  use  in  the  beamformer.  It 
calculates  60  second  averages  for  all  data  fields,  converts 
pressure  to  depth  and  calculates  tilt  direction  based  on 
current  velocity. 


/wMficit 

*  PROGRAM  SACM1.C  vsn  1.0 

*  WRITTEN  BY:  Steven  Crocker 

*  LAST  UPDATE:  August  6,  1991 
* 

*  This  program  takes  SACM  data  from  the  Heard  Island  West  Coast 

*  Array  (upper  instrument  package)  and  condenses  it. 
* 

*  Pressure  is  converted  to  depth. 
* 

*  The  conductivity  is  not  processed.  All  values  output  are  60 

*  second  averages  of  the  input. 

* 

A  A  A  A  A  / 


fldefine  VERSION  ANSI 
^define  SELECT  OUTPUT  OFF 
^define  C  25 
#define  inBUFFER  1200 
tfdefine  outBUFFER  50 
#define  PI  3.14159265359 
#define  RADIAN  57.2957795131 

#include<stdio.h> 


/*  either  ANSI  or  UNIX  */ 

/*  ON  or  OFF  */ 

/*  max  length  for  file  names  */ 

/*  input  buffer  size  */ 

/*  output  buffer  size  */ 
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^include  <malloc.h> 
#include<math.h> 

#lf  (VERSION==ANSI) 

tt  i  nc  lude<std  Lib.  h> 

void  ExitOnErrorCchar  error_txt[]); 

tfendif 

mainO 
i 

FILE  *sacm_fp,  *out_fp; 

char  c,  sacm_f  i  leCC],  out_f  i  leCG; 

int  i,  j,  count,  loop_count,  inBUFFER_fulL,  out_f LagC6]; 

float  t1,  avgTilt,  avgPress,  avgVN,  avgVE,  avgTemp,  avgAngle, 
avgVelocity,  junk; 

static  float  inTiltCinBUFFER],  inPressCinBUFFER],  inVNCinBUFFER], 
inVECinBUFFER],  inTimeCinBUFFER],  inTempCinBUFFER], 
outTimeCoutBUFFER],  outTi ItCoutBUFFER], 
outPressCoutBUFFERD,  outAngleLoutBUFFER], 
outVelocityEoutBUFFER],  out TempCout BUFFER] ; 

printf ("\n\n\nEnter  file  name  for  SACM  data.Xn"); 

scanf ("%s",sacm_fi le); 

printf ("\n\n"); 

printf ("Enter  file  name  for  output. \n"); 

scanf ("%s"/out_fi le); 

printf ("\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n"); 

printf ("The  following  pertains  to  output  file  format .\n\n"); 

printf ("Units  of  measure:   Time  index  is  minutes. \n"); 
printfC  Angles  in  degrees\n"); 

printfC  Compass  direction  (TRUE)\n"); 

printfC  Velocity  in  meters/second\n"); 

printfC  Temperature  in  degrees  C\n"); 

printfC  Depth  in  meters\n\n"); 

printf ("NOTE:   Output  column  order  is  same  as  above. \n\n\n"); 

for  (i=0;i<6;i++)  out_f lag[i]=1; 

#if  (SELECT_OUTPUT==ON) 

printf ("Include  TIME?  (y  or  n)\n"); 

scanf ("%s",&c); 

if(c==89  !!  c==121)  out_flag[0]=1; 

else  out_f lag[0]=0; 

printf("\n\n"); 

printf ("Include  TILT?  (y  or  n)\n"); 

scanf ("%s",&c); 

if(c==89  Jl  c==121)  out_flag[1]=1; 

else  out_f lag[1]=0; 

printf("\n\n"); 

printfCInclude  CURRENT  DIRECTION?  (y  or  n)\n"); 

scanf ("%s",&c ); 

if(c==89  Jl  c==121)  out_f lagC2]=1; 

else  out_f lag[2]=0; 

printf("\n\n"); 

printfCInclude  CURRENT  VELOCITY?  (y  or  n)\n"); 

scanf ("%s" ,&c); 

if(c==89   ;|    c==121)  out_flag[3>1; 

else  out_flag[33=0; 

printf ("\n\n"); 

printfCInclude  TEMPERATURE?   (y  or  n)\n"); 
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scanf("%s",&c); 

if(c==89   ||    c==121)  out_fLag[A>1; 

else  out_f  lag [411=0; 

printf ("\n\n"); 

printfC'IncLude  DEPTH?   (y  or  n)\n"); 

scanf("%s",&c); 

if(c==89   ||   c==121)  out_flag[5>1; 

else  out_f LagC5]=0; 

tfendif 

inBUFFER_full=inBUFFER; 
loop_count=0; 

if   ((sacm_fp=fopen(sacm_file,    "rt"))    !=  NULL)   ; 

else  ExitOnErrorC'Error  opening  input  data  file"); 

if  ((out_fp=fopen(out_file,  "wt"))  !=  NULL)  ; 

else  ExitOnErrorC'Error  opening  output  data  file"); 

while  (1) 
i 

for  (i=0;i<inBUFFER;i++) 
{ 

if  (fscanf(sacm_fp/  "%f  %f  %f  7.1   %f  %f  %f\n", 

&inTime[i],&inVN[i],  SinVEEiD,  SinTempEi], 
&inTilt[i],  SinPressCi],  &junk)  !=  EOF)  ; 

else  inBUFFER_full=i-2; 

>  /*  for  */ 
j=0; 

for  (i=0;i<inBUFFER_full;i++) 
i 

t1=inTime[i]; 

avgTilt=0.0; 

avgPress=0.0; 

avgVN=0.0; 

avgVE=0.0; 

avgTempO.O; 

count=0; 

while((inTimeCi]<t1+60.0)  &&  Ci<inBUFFER  full)) 

< 

avgTilt=avgTilt  +  inTi  ltCi]/10.0; 

avgPress=avgPress  +  inPressli]*1000.0; 

avgVN=avgVN  +  i nVN[i]/1 000.0; 

avgVE=avgVE  +  inVE[i]/1000.0; 

avgTemp=avgTemp  +  inTemp[i 3/100.0; 

i++; 

count++; 
>  /*  while  */ 

avgTi lt=avgTi lt/(f loat)count; 
avgPress=avgPress/(f  loat)count; 
avgTemp=avgTemp/ ( f loat ) count ; 
avgVN=avgVN/ ( f loat ) count; 
avgVE=avgVE/(f loat) count ; 
avgVelocity=sqrt(avgVN*avgVN+avgVE*avgVE); 
avgAngle=atan2(avgVE,   avgVN)*RADIAN; 
if   (avgAngle<=0.0)  avgAngle=360.0+avgAngle; 
outTime[j]=(f  loat)(  j+D+outBUFFERnoop^oun^- 
outTi  lt[j>avgTi It; 

outPressC jXavgPress-101 325. 0)/(9. 80665*1026.0); 
outTempE j ]=avgTemp; 
outVe  loci ty[j]=avgVe loci ty ; 
out Ang leC j  3=avgAng I e; 

j++; 

>  /*  for  */ 

if  (inBUFFER_full  !=  inBUrFER)  j— ; 

for  (i=0;i<j;i++) 

{ 

if  (out_f lag[0])  fprintf (out_fp,  "%6.0f",  outTimeEi]), 
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if  (out_f lag[1])  fprintf (outjfp, 

if  (out_f LagC2])  fprintf (outjfp, 

if  (out_f LagC3])  fprintf (out_fp, 

if  (out_f lag[4])  fprintf (outjfp, 

if  (outjf  lag[53)  fprintf (out_fp, 
fprintf(out_fp,  "\n"); 
/*  for  */ 


if  (inBUFFER_fuLL  != 

Loop_count++; 
>  /*  while  */ 
exitd); 


inBUFFER)  exit(O); 


X8.5f",  outTiltCi]); 
%8.3f ",  outAngleCiD); 
%8.6f",  outVeLocity[i]); 
%8.5f",  outTempCi]); 
X8.0f",  outPressCi]); 


#if  (VERSION==ANSI) 

void  ExitOnErrorCchar  error_txt[]) 

#elif  (VERSION==UNIX) 

ExitOnError(error_txt) 

char  error_txt[3; 

#endif 

{. 

fprintf (stderr, "Program  run-time  error  ...\n"); 

fprintf (stderr, "%s\n",error_txt); 

fprintf (stderr, ".. .now  exiting  to  system. . .\n"); 

exitd); 


D .       SACM2 . C 


This  program  takes  the  output  of  sacral. c  as  input.  It 
locates  and  copies  a  user  defined  subset  of  the  tilt  data  for 
use  in  beamforming. 


/jtiricick 

*      PROGRAM  SACM2.C  vsn  1.0 
WRITTEN  BY:  Steven  Crocker 
LAST  UPDATE:  August  6,    1991 


This  program  takes  data  processed  by  SACM1 . C  and  cuts 
a  user  defined  segment  from  it.  The  segment  is  retained 
as  the  output  data  file.  The  input  file  is  not  affected. 
The  output  time  base  may  either  be  normalized  to  1  or  may 
retain  the  original  values. 

The  option  is  given  to  accept  a  default  output  format.  This 
output  format  coincides  with  the  required  input  format  to  the 
time  domain  modal  beamformer. 

The  output  elements  are  selectable  if  desired. 
/ 


/*  either  ANSI  or  UNIX      */ 
/*  max  length  for  file  names  */ 


^define  VERSION  ANSI 
fldefine  C  25 
#include<stdio.h> 
^include  <malloc.h> 


#if    (VERSION==ANSI) 

#include<stdlib.h> 

void   ExitOnErrorCchar  error_txt[]); 

flendif 
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mainO 
i 


FILE  *input_fp,  *out_fp; 

char  c,  Default,  input_f i LeCC],  out_f i le[C]; 

int  i=0,  outl,  outN,  out_f lag[7]; 

float  Tilt,  Depth,  Time,  Temp,  Angle,  Velocity; 

printf ("\n\n\nEnter  file  name  for  input  data.Nn"); 
scanf ("%s",input_f i le); 

if  ((input_fp=fopen(input_file,  "rt"))  !=  NULL)  ; 
else  ExitOnErrorC'Error  opening  input  file"); 
printf("\n\n"); 

printf ("Enter  file  name  for  output. \n"); 
scanf ("%s",out_fi le); 

if  ((out_fp=fopen(out_file,  "wt"))  !=  NULL)  ; 

else  ExitOnErrorC'Error  opening  output  file"); 

printf ("\n\n\nDo  you  want  the  default  output  conf iguration?\n"); 

printf ("This  option  formats  the  output  for  use  in  the  beamformer.Nn"); 

scanf ( "%s" , Def au It ) ; 

printf (" \n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n"); 

if(Default  !=  'n'  S&  Default  !=  'N') 

{. 

printf ("The  following  pertains  to  output  file  format. \n\n"); 

printf ("Units  of  measure:   Time  Index  is  minutes. \n"); 

pnntfC  Angles  in  degrees\n"); 

printfC  Compass  direction  (TRUE)\n"); 

printfC  Velocity  in  meters/second\n"); 

printfC  Temperature  in  degrees  C\n"); 

printfC  Depth  in  meters\n\n"); 

printf ("NOTE:  Output  column  order  is  same  as  above. \n\n\n"); 

printf ("Include  TIME?  (y  or  n)\n"); 
scanf  V"/>s",c); 

if(c==89  |!  c==121)  out_flag[0]=1; 
else  out_f lag[0D=0; 
printf ("\n\n"); 

printfC'Include  TILT?  (y  or  n)\n"); 

scanf  C7.s",c  ); 

if(c==89  |l  c==121)  out_flag[1]=1; 

else  out_f lag[1]=0; 

printf ("\n\n"); 

printfC'Include  CURRENT  DIRECTION?  (y  or  n)\n"); 

scanf ("%s",c ); 

if(c==89  ||  c==121)  out_flag[2]=1; 

else  out_f lag[2]=0; 

printf ("\n\n"); 

printfC'Include  CURRENT  SPEED?  (y  or  n)\n"); 

scanf ("%s",c); 

if(c==89  ||  c==121)  out_flag[3]=1; 

else  out_f lag[3]=0; 

printf("\n\n"); 

printfC'Include  TEMPERATURE?  (y  or  n)\n"); 

scanf ("%s",c); 

if(c==89  ||  c==121)  out_flag[4]=1; 

else  out_f lagCA]=0; 

printf("\n\n"); 

printfC'Include  DEPTH?  (y  or  n)\n"); 
scanf ("%s",c); 

if(c==89  ||  c==121)  out_flagC5]=1; 
else  out_f lag[5]=0; 
>  /*  if  */ 
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else 
{ 

outjf  lag[0>0; 
out_f lag[1]=1; 
out_flag[2]=1; 
out_fLagC3]=0; 
out_f lag[4]=0; 
out_f LagC5]=1; 

>  /*  else  */ 

printf ("\n\nWhat  time  index  of  the  input  file  should  be  the  first\n"); 

printf ("element  of  the  output  file?\n"); 

scanf("7.i",&out1); 

printf("\n\n"); 

printf ("How  many  elements  should  the  output  file  contain?\n"); 

scanf("%i",&outN); 

printf("\n\n"); 

if  (Default  !=  'n'  &&  Default  !=  'N') 
i 

printf ("Should  the  time  base  be  normalized  to  begin  at  t=1?\n"); 

scanf ("%s",c); 

if(c==89  il  c==121)  out_flagC6]=1; 

else  out_f lag[6]=0; 

printf ("\n\n\n"); 

>  /*  if  */ 


while  (1) 
{ 


if  (fscanf(input_fp,  "%f  %f  %f  %f  %f  %f  \n", 

STime,  &Tilt,  SAngle,  SVelocity,  STemp,  SDepth)  !=  EOF) 

{ 

if((Time  >=   (float)outD  &&   (Time  <  (f Loat)(out1+outN))) 
{ 

if(out_flag[0]  &S   !out_f lag[6]) 
fprintf (out_fp,"%f ",   Time); 

else  if   (out_f lagCO]  &&  out_flag[6]) 

fprintf(out_fp,    "%f",    (f loat)(i+D); 


if  (out_f lagCID)  fprintf (out_fp,  " 

%f 

'/  Tilt); 

if  (out_f lagC2])  fprintf (out_fp,  " 

'/.1 

',  Angle); 

if  (out_f lagC3])  fprintf (out_fp,  " 

%f 

',  Velocity); 

if  (out_f lag[4])  fprintf (out_fp,  " 

%f 

',  Temp); 

if  (out_f lag[5])  fprintf (outjfp,  " 

%f 

',  Depth); 

fprintf(out  fp,  "\n"); 

i++  ; 

}  /*  if  */ 

else  if (Time  >=  (f loat)(out1+outN)) 
{ 

fclose(input_fp); 

fclose(out_fp); 

exit(O); 

>  /*  else  if  */ 

>  /* 

if  */ 

else 
{ 

f close(input_fp); 
f close(out_fp); 

ExitOnErrorC'EOF  encountered  in  input  data 

fil 

e"); 

y  /* 

else  */ 

>  /*  while 

*/ 

exitd); 

#if    (VERSION==ANSI) 
void   ExitOnError(char   error_txt[]) 
#elif    (VERSION==UNIX) 
ExitOnError(error   txt) 
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char  error_txt[]; 

#endif 

{ 

fprintf (stderr, "Program  run-time  error  ...\n"); 

fprintf (stderr, "%s\n",error_txt); 

fprintf (stderr,". . .now  exiting  to  system. . .\n"); 

exit(l); 
> 

E.   ARRAYTEST.C 

This  program  performs  various  statistical  tests  on  each 
element  of  the  acoustic  array.  It  is  useful  for  isolating 
hydrophones  which  have  failed  during  the  experiment.  The 
statistical  tests  are  an  adaptation  of  those  found  in 
Reference  11. 


/'k'k'kic 

*  PROGRAM:  ARRAYTST  vsn  1.0 

*  WRITTEN  BY:  Steven  Crocker 

*  LAST  UPDATE:  October  21,  1991 


I 

#def ine 
#def ine 
#def ine 


UNIX  VERSION 
CHANNELS  32 
F  SAMPLE  228 


/*  either  ANSI  or  UNIX  */ 

/*  number  of  hydrophones  on  array  */ 

/*  sampling  frequency  */ 

/*  duration  of  time  averages  */ 


#define  BUFFER_TIME  60 

#include<stdio.h> 
#i  nc lude<ma  I loc .  h> 
#include<math.h> 


#if  defined  (  ANSI  ) 

#include<f loat.h> 

#include<stdlib.h> 

int  getlnput(void); 

int  momentsCf Loat  *data,  int  n,  float  *ave,  float  *ave2,  float  *adev, 

float  *sdev,  float  *svar,  float  *skew,  float  *curt) 
float  *vector(int  length); 
float  **matrix(int  row,  int  col); 
int  f ree_matrix(f loat  **m,  int  row); 
void  ExitOnErrorCchar  error_txt[]); 
#elif  defined  (  UNIX  ) 
int   vectorO, 

matrixO, 

getlnputO, 

momentsO, 

f ree_matrix(), 

ExitOnErrorO; 
#endif 

/*  Global  Variable  Declarations  */ 
float  **inS0UND; 
int  f irstBuffer=1; 

mainO 

int  i,  j=1,  k,  flagC8]; 
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float  ave,  ave2,  adev,  sdev,  svar,  skew,  curt,  *temp; 

char  fileNameCM)],  key; 

FILE  *fp1,  *fp2,  *fp3,  *fp4,  *fp5,  *fp6,  *fp7; 

inSOUND=( float  **)matrix(CHANNELS/  BUFFER_TIME*F_SAMPLE); 
temp=(float  *)vector(BUFFER_TIME*F_SAMPLE); 

for  (i=0;i<=7;i++)  flag[i]=0; 

printf ("Answer  y  or  n  to  the  following  output  options. \n"); 

printf ("Average  value  (  over  %i  seconds  ):  ",  BUFFER_TIME); 

scanf("%s",  ftkey); 

if(key==89  ||  key==121) 

{ 

flagC1]=1; 

printf ("Enter  file  name:  "); 

scanf("%s",  fileName); 

printf ("\n\n"); 

if  ((fp1=fopen(fileName,  "wt"))  ==  NULL) 
ExitOnErrorC'Error  opening  file."); 

>  /*  if  */ 

printf ("Mean  squared  value  (  over  %i  seconds  ):  ",  BUFFER_TIME); 
scanf("%s",  &key); 
if(key==89  ||  key==121) 
•C 

flagC2]=1; 

printf ("Enter  file  name:    "); 

scanfC/is",   fileName); 

printf("\n\n"); 

if  ((fp2=fopen(fileName/  "wt"))  ==  NULL) 
ExitOnErrorC'Error  opening  file."); 

>  /*  if  */ 

printf ("Average  deviation  (  over  %i  seconds  ):  ",  BUFFER_TIME); 
scanf("%s",  &key); 
if(key==89  j |  key==121) 
{ 

flag[3]=1; 

printf ("Enter  file  name:  "); 

scanf("%s",  fileName); 

printf ("\n\n"); 

if  ((fp3=fopen(fileName,  "wt"))  ==  NULL) 
ExitOnErrorC'Error  opening  file."); 

>  /*  if  */ 

printf ("Standard  deviation  (  over  %i  seconds  ):  ",  BUFFER_TIME); 
scanfC7.s",  &key); 
if(key==89  | |  key==121) 
{ 

f  lag[A>1; 

printf ("Enter  file  name:  "); 

scanf("%s",  fileName); 

printf("\n\n"); 

if  ((fp4=fopen(fileName,  "wt"))  ==  NULL) 
ExitOnErrorC'Error  opening  file."); 

>  /*  if  */ 

printf ("Variance  (  over  %i  seconds  ):  ",  BUFFER_TIME); 
scanf('7.s",  &key); 
if(key==89  ||  key==121) 
{ 

flagC5]=1; 

printf ("Enter  file  name:  "); 

scanf("%s",  fileName); 

printf("\n\n"); 

if  ((fp5=fopen(fileName,  "wt"))  ==  NULL) 
ExitOnErrorC'Error  opening  file."); 

>  /*  if  */ 
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printfC'Skewness  (  over  %i  seconds  ):  ",  BUFFER_TIME); 
scanf("%s",  &key); 
if(key==89  ||  key==121) 
i 

flagC6]=1; 

printf ("Enter  file  name:  "); 

scanf("%s",  fileName); 

printf("\n\n"); 

if  ((fp6=fopen(fileName,  "wt"))  ==  NULL) 
ExitOnErrorC'Error  opening  file."); 
>  /*  if  */ 


printf ("Kurtosis  (  over  %i  seconds  ):  ",  BUFFER_TIME); 
scanf("%s",  Skey); 
if(key==89  ||  key==121) 
{ 

fLag[7]=1; 

printf ("Enter  file  name:  "); 

scanf("%s",  fileName); 

printf("\n\n"); 

if  ((fp7=fopen(fileName,  "wt"))  ==  NULL) 
ExitOnErrorC'Error  opening  file."); 
>  /*  if  */ 


while  (j<120) 

i 

get  Input (); 

if  (flagCI])  fprintf(fp1,"%i 

if  (flag[2])  fprintf (fp2,"%i 

if  (flag[3])  fprintf (fp3,'7.i 

if  (flag[4])  fprintf (fp4,"%i 

if  (flag[53)  fprintf (fp5,"%i 

if  (flagC63)  fprintf  (fp6,"°/i 

if  (flag[7])  fprintf (fp7,"%i 


",j>; 
\j); 
",j>; 
",j>; 
",j>; 
",j>; 
",j); 


for  (i=1;i<=CHANNELS;i++) 
i 

for(k=1;k<=BUFFER_TIME*F_SAMPLE;k++) 
temp[k:=inSOUNDCi]Ck]; 

moments(temp/BUFFER_TIME*F_SAMPLE,&ave,&ave2, 

Sadev^sdev^svar/Sskew^curt); 
if  (flag[1])  fprintf(fp1, 
if  (flagC23)  fprintf(fp2, 
if  (flag[33)  fprintf(fp3, 
if  (flag[4])  fprintf(fp4, 
if  (flag[51)  fprintf(fp5, 
if  (flag[61)  fprintf(fp6, 
if  (flag[7])  fprintf(fp7, 
>  /*  for  */ 


1  %f  ' 

,   ave); 

1  %f  ' 

,   ave2); 

1  %f  ' 

',  adev); 

1  %f  ' 

',   sdev); 

'  %f 

',    svar); 

"  %f 

',    skew); 

"  %f 

',    curt); 

if  (flag[1])  fprintf(fp1,"\n"); 

if  (flagC2])  fprintf (fp2,"\n"); 

if  (flag[3])  fprintf (fp3,"\n"); 

if  (flag[43)  fprintf Cfp4,"\n"); 

if  (flag[53)  fprintf (fp5,"\n"); 

if  (flag[6])  fprintf (fp6,"\n"); 

if  (flag[7])  fprintf (fp7,"\nM); 
firstBuffer=0; 

>  /*  while  */ 
exit(O); 

>/* 


*   FUNCTION:  getlnputO 
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*  This  function  handles  all  acoustic  input.   It  also  provides  one  of 

*  two  normal  process  terminations  available  in  the  program. 
* 


ExitOnErrorO 


*  Arguments: 
* 

none 

*  Return  value: 
* 

*  Functions  called: 

0 

vectorO 

*  Definitions  called: 

ANSI 

* 

F  SAMPLE 

* 
* 
*  Global  variables  called: 

BUFFERJTIME 

inSOUND[][] 

*  Significant  memory  allocation: 

* 

A  A  A  A  A  / 

diskBufferC] 

#if  defined  (  ANSI  ) 

int  getlnput(void) 

#elif  defined  (  UNIX  ) 

get  Input () 

flendif 

{ 

int  i,  j,  items,  buffer; 

float  *diskBuffer; 

char  fileName[803; 

static  FILE  *fpStatic; 

UNIX 
CHANNELS 


f irstBuffer 


if  (f irstBuffer) 
{ 

printf ("Enter  file  name  for  input  acoustic  data:  "); 

ifUscanf  ("%s",   f  i  leName))==EOF) 

ExitOnErrorC'Fatal  error   in  scanfO"); 

printf ("\n\n"); 

if  ((fpStatic=fopen(fileName/  "rb"))  ==  NULL) 

ExitOnErrorC'Error  opening  INPUT  ACOUSTIC  data  file."); 

>  /*  if  */ 

/*  Memory  Allocation  */ 

buffer=BUFFER_TIME*F_SAMPLE*CHANNELS; 
diskBuffer=(f loat*) vector (buffer ); 

items=fread((char*)(diskBuffer+1)/  si zeof (float),  buffer,  fpStatic); 
if  (items==buffer)  /*continue*7; 
else  if (ferror(fpStatic)  !=  0) 

ExitOnErrorC'Error  encountered  while  reading  input  acoustic  data"); 
else  if (feof (fpStatic)  !=  0) 
i  ,____ 

printf ("\t    End  of  File  reached:  EXECUTION  COMPLETED"); 

f close(fpStatic); 
exit(0); 

>  /*  else  if  */ 

else  ExitOnErrorC'Unknown  error  handling  acoustic  input  file."); 

for  (i=1;  i<=BUFFER_TIME*F_SAMPLE;  i++) 
{ 

for  (j=1;  j<=CHANNELS;  j++) 

inS0UND[j][i]  =  disk3uff er[CHANNELS*(i-1 )+j3; 
>  /*  for  */ 

/*  Deallocate  Memory  */ 
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f ree((char*)diskBuffer); 
returnC  0  ); 

/*****  END  get  Input  *****/ 

/a***********************/ 

/ 


*  FUNCTION:  moments () 
* 

*  This  function  calculates  the  mean,  mean  squared  value,  average 

*  deviation,  standard  deviation,  variance,  skewness  and  kurtosis 

*  of  a  data  vector. 


data 

ave 

adev 

svar 

curt 


*  Arguments: 
* 
* 
* 
* 
* 

*  Return  value: 
* 

*  Functions  called: 
* 

*  Definitions  called: 
* 

*  Global  variables  called: 
* 

*  Significant  memory  allocation:   none 
* 


n 

ave2 
sdev 
skew 


ExitOnErrorO 

ANSI 

none 


UNIX 


#if  defined  (  ANSI  ) 

int  moments (float  *data,  int  n,  float  *ave,  float  *ave2,  float  *adev, 

float  *sdev,  float  *svar,  float  *skew,  float  *curt) 
flelif  defined  (  UNIX  ) 

moments (data, n, ave, ave2, adev, sdev, svar, skew,  curt) 
int  n; 

float  *data,  *ave,  *ave2,  *adev,  *sdev,  *svar,  *skew,  *curt; 
#endif 
{ 

int  j; 

float  s,p; 

if  (n<=1)  ExitOnErrorC'n  must  be  at  least  2  in  moment ()"); 

s=0.0; 

*ave2=0.0; 

for  (j=1;j<=n;j++) 

i 

s  +=  dataEj]; 

*ave2  +=  dataCj]*data[j3; 

>  /*  for  */ 
*ave=s/n; 
*ave2  /=  n; 

*adev=(*svar)=(*skew)=(*curt)=0.O; 
for  (j=1;j<=n;j++) 

{ 

*adev  +=  fabs(s=data[j]-(*ave)); 
*svar  +=  (p=s*s); 
*skew  +=  (p  *=  s); 
*curt  +=  (p  *=  s); 

>  /*  for  */ 
*adev  /=  n; 
*svar  /=  (n-1); 
*sdev=sqrt (*svar ) ; 
if  (*svar) 

i 

•skew  /=  (n*(*svar)*(*sdev)); 
*curt=(*curt)/(n*(*svar)*(*svar))-3.0; 

>  /*  if  */ 

else  ExitOnErrorC'No  skew/kurtosis  when  variance  =  0  (in  momentO)"), 
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return(  0  ); 
/*****  END  moment  aaaaa/ 

/A  A  A  A  A 

*  FUNCTION:  vectorO 
* 

*  This  function  allocates  memory  for  UNIT  OFFSET  vectors. 

* 


length 

*v 

ExitOnErrorO 

ANSI 

none 


*  Arguments: 
* 

*  Return  value: 
* 

*  Functions  called: 
* 

*  Definitions  called:  ANSI  UNIX 
* 

*  Global  variables  called: 
* 

*  Significant  memory  allocation:   v[] 
* 

#if  defined  (  ANSI  ) 

float  *vector(int  length) 

#elif  defined  (  UNIX  ) 

vectorC length) 

int  length; 

flendif 

{ 

float  *v; 

if  ((v=(f loat*)ma Hoc (dength+1  )*si zeof (float )))==NULL) 
ExitOnErrorC'Memory  allocation  failure  in  vectorO."); 
#if  defined  (  ANSI  ) 

return  v; 
#elif  defined  (  UNIX  ) 

return  (long  int)v; 
#endif 

y    /irkirft A A  A"  A A AAAAAAAAA"AAAA/ 

/*•***  END  vector  •*•**/ 

<  i    i    i    i    ■    ■    i    i    i    i    i    ■    t    ■    i    i    ■    ■    i    i    ■    i     I 

firkickit 

*  FUNCTION:  matrixO 
* 

*  This  function  allocates  memory  for  UNIT  OFFSET  2-D  arrays. 
* 

col 


UNIX 


*  Arguments: 
* 

row 

*  Return  value: 

**m 

*  Functions  called: 

ft 

ExitOnErrorO 

*  Definitions  called: 

ANSI 

*  Global  variables  called: 

none 

* 

*  Significant  memory  allocation: 

ft 

m[][] 

'k'k'k'k'k  / 

#if  defined  (  ANSI  ) 

float  **matrix(int  row,  int  col) 

#elif  defined  (  UNIX  ) 

matrix(row,col) 

int  row,  col; 

flendif 

i 

int  i; 

float  **m; 
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if  ((m=(float**)ma I loc( (unsigned) (row+1)*sizeof(float*)))==NULL) 

ExitOnErrorC'Allocation  failure  1  in  matrixO"); 
for  (i=1;  i<=row;  i++) 
i 

if  ((m[i]=(f loat*)ma Hoc ((unsigned) (col+1 )*sizeof (float )))==NULL) 
ExitOnErrorC'Allocation  failure  2  in  matrixO"); 
>  /*  for  */ 
#if  defined  (  ANSI  ) 

return  m; 
tfelif  defined  (  UNIX  ) 
return  (long  int)m; 
#endif 

>  /AAAAAAAAAAAAAAAAAAAAAA/ 
/*****  END  matrix  *****/ 
/AAAAAAAAAAAAAAAAAAAAAA/ 
/AAA  A  A 

*  FUNCTION:  free_matrix() 
* 

*  This  function  deallocates  memory  from  UNIT  OFFSET  2-D  arrays. 
* 

*  Arguments:  mCDC]  row 
* 

*  Return  value:  0 
* 

*  Functions  called:  none 

A 

A  Definitions  called:  ANSI  UNIX 

A 

a  Global  variables  called:        none 

A 

*  Significant  memory  allocation:   none 

A 


/ 

#if  defined  (  ANSI  ) 

void  f ree_matnx(f loat  **m,  int  row) 

#elif  defined  (  UNIX  ) 

f ree_matrix(m,  row) 

float  **m; 

i  nt  row; 

flendif 

{ 

int  i; 


for(i=row;  i>=1;  i — ) 
f ree((charA)m[i]); 
f ree((charA)m); 
return(  0  ); 

/a**aa  END  free  matrix  a*a*a/ 
/irkirkirkirkiciciciclrkiriciclcirjricicirkirirk  j 
/icirkirk 

*  FUNCTION:  Exi tOnErrorO 

A 

A  This  function  performs  an  abnormal  process  termination. 
A 

*  Arguments:  error_txt[] 

A 

*  Return  value:  none 

A 

*  Functions  called:  none 
A 

A  Definitions  called:  ANSI  UNIX 

A 

*  Global  variables  called:        none 

A 

*  Significant  memory  allocation:   none 

A 

Irk/ 
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#if  defined  (  ANSI  ) 

void  ExitOnErrorCchar  error_txt[]) 

#elif  defined  (  UNIX  ) 

ExitOnError(error_txt) 

char  error_txt[]; 

#endif 

fprintf (stderr, "Program  run-time  error  .. 
fprintf (stderr, "%s\n",error_txt); 
fprintf (stderr,". . .now  exiting  to  system. 
exit(O); 

/*****  END  ExitOnError  *****/ 


An"); 
An"); 
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